Metrics for assessing community assembly processes
source('../env.R')community_data = read_csv(filename(COMMUNITY_OUTPUT_DIR, 'community_assembly_metrics_using_relative_abundance.csv'))Rows: 336 Columns: 43── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
dbl (43): mntd_normalised, mntd_actual, mntd_min, mntd_max, mntd_mean, mntd_sd, fdiv_normalised, fdiv_actual, fdiv_min, fdiv_max, fdiv_mean, fdiv_sd, mass_fdi...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.head(community_data)colnames(community_data) [1] "mntd_normalised" "mntd_actual" "mntd_min" "mntd_max"
[5] "mntd_mean" "mntd_sd" "fdiv_normalised" "fdiv_actual"
[9] "fdiv_min" "fdiv_max" "fdiv_mean" "fdiv_sd"
[13] "mass_fdiv_normalised" "mass_fdiv_actual" "mass_fdiv_min" "mass_fdiv_max"
[17] "mass_fdiv_mean" "mass_fdiv_sd" "locomotory_trait_fdiv_normalised" "locomotory_trait_fdiv_actual"
[21] "locomotory_trait_fdiv_min" "locomotory_trait_fdiv_max" "locomotory_trait_fdiv_mean" "locomotory_trait_fdiv_sd"
[25] "trophic_trait_fdiv_normalised" "trophic_trait_fdiv_actual" "trophic_trait_fdiv_min" "trophic_trait_fdiv_max"
[29] "trophic_trait_fdiv_mean" "trophic_trait_fdiv_sd" "gape_width_fdiv_normalised" "gape_width_fdiv_actual"
[33] "gape_width_fdiv_min" "gape_width_fdiv_max" "gape_width_fdiv_mean" "gape_width_fdiv_sd"
[37] "handwing_index_fdiv_normalised" "handwing_index_fdiv_actual" "handwing_index_fdiv_min" "handwing_index_fdiv_max"
[41] "handwing_index_fdiv_mean" "handwing_index_fdiv_sd" "city_id" Join on realms
city_to_realm = read_csv(filename(CITY_DATA_OUTPUT_DIR, 'realms.csv'))Rows: 337 Columns: 2── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): core_realm
dbl (1): city_id
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.community_data_with_realm = left_join(community_data, city_to_realm)Joining with `by = join_by(city_id)`Cities as points
city_points = st_centroid(read_sf(filename(CITY_DATA_OUTPUT_DIR, 'city_selection.shp'))) %>% left_join(community_data_with_realm)Warning: st_centroid assumes attributes are constant over geometriesWarning: st_centroid does not give correct centroids for longitude/latitude dataJoining with `by = join_by(city_id)`city_points_coords = st_coordinates(city_points)
city_points$latitude = city_points_coords[,1]
city_points$longitude = city_points_coords[,2]world_map = read_country_boundaries()Reading layer `WB_countries_Admin0_10m' from data source `/Users/james/Dropbox/PhD/WorldBank_countries_Admin0_10m/WB_countries_Admin0_10m.shp' using driver `ESRI Shapefile'
Simple feature collection with 251 features and 52 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -180 ymin: -59.47275 xmax: 180 ymax: 83.6341
Geodetic CRS: WGS 84Warning: st_simplify does not correctly simplify longitude/latitude data, dTolerance needs to be in decimal degreesLoad community data, and create long format version
communities = read_csv(filename(COMMUNITY_OUTPUT_DIR, 'communities_for_analysis.csv'))Rows: 2427 Columns: 12── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): city_name, jetz_species_name, seasonal, presence, origin
dbl (4): city_id, distance_to_northern_edge_km, distance_to_southern_edge_km, relative_abundance_proxy
lgl (3): present_urban_high, present_urban_med, present_urban_low
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.communitiescommunity_summary = communities %>% group_by(city_id) %>% summarise(regional_pool_size = n(), urban_pool_size = sum(relative_abundance_proxy > 0))
community_summaryLoad trait data
traits = read_csv(filename(TAXONOMY_OUTPUT_DIR, 'traits_jetz.csv'))Rows: 304 Columns: 6── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): jetz_species_name
dbl (5): gape_width, trophic_trait, locomotory_trait, mass, handwing_index
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.head(traits)Load realm geo
resolve = read_resolve()
head(resolve)Simple feature collection with 6 features and 15 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -162.1547 ymin: -69.55876 xmax: 158.6167 ymax: 61.53428
Geodetic CRS: WGS 84Summary metrics by Realm
test_required_values = function(name, df) {
cat(paste(
test_value_wilcox(paste(name, 'MNTD'), df$mntd_normalised),
test_value_wilcox(paste(name, 'Beak Gape FDiv'), df$gape_width_fdiv_normalised),
test_value_wilcox(paste(name, 'HWI FDiv'), df$handwing_index_fdiv_normalised),
test_value_wilcox(paste(name, 'Mass FDiv'), df$mass_fdiv_normalised),
nrow(df),
sep = "\n"))
}test_required_values('Global', community_data_with_realm)Global MNTD median 0.49
Global Beak Gape FDiv median 0.59 ***
Global HWI FDiv median 0.66 ***
Global Mass FDiv median 0.65 ***
336unique(community_data_with_realm$core_realm)[1] "Oceania" "Nearctic" "Neotropic" "Palearctic" "Afrotropic" "Indomalayan" "Australasia"test_required_values('Nearctic', community_data_with_realm[community_data_with_realm$core_realm == 'Nearctic',])Warning: cannot compute exact p-value with tiesWarning: cannot compute exact p-value with tiesNearctic MNTD median 0.65
Nearctic Beak Gape FDiv median 0.63
Nearctic HWI FDiv median 0.23 **
Nearctic Mass FDiv median 0.48
66test_required_values('Neotropic', community_data_with_realm[community_data_with_realm$core_realm == 'Neotropic',])Neotropic MNTD median 0.5
Neotropic Beak Gape FDiv median 0.47
Neotropic HWI FDiv median 0.52
Neotropic Mass FDiv median 0.66 ***
64test_required_values('Palearctic', community_data_with_realm[community_data_with_realm$core_realm == 'Palearctic',])Palearctic MNTD median 0.59 *
Palearctic Beak Gape FDiv median 0.93 ***
Palearctic HWI FDiv median 0.4
Palearctic Mass FDiv median 0.55
74test_required_values('Afrotropic', community_data_with_realm[community_data_with_realm$core_realm == 'Afrotropic',])Afrotropic MNTD median 0.16 *
Afrotropic Beak Gape FDiv median 0.42
Afrotropic HWI FDiv median 0.54
Afrotropic Mass FDiv median 0.34
9test_required_values('Indomalayan', community_data_with_realm[community_data_with_realm$core_realm == 'Indomalayan',])Indomalayan MNTD median 0.44 ***
Indomalayan Beak Gape FDiv median 0.46
Indomalayan HWI FDiv median 0.88 ***
Indomalayan Mass FDiv median 0.8 ***
116test_required_values('Australasia', community_data_with_realm[community_data_with_realm$core_realm == 'Australasia',])Australasia MNTD median 0.53
Australasia Beak Gape FDiv median 0.46
Australasia HWI FDiv median 0.78
Australasia Mass FDiv median 0.55
6What families exist in which realms?
communities %>%
left_join(city_to_realm) %>%
mutate(family = gsub( "_.*$", "", jetz_species_name)) %>%
dplyr::select(family, core_realm) %>%
distinct() %>%
arrange(core_realm)Joining with `by = join_by(city_id)`Summary metrics by introduced species
communities = read_csv(filename(COMMUNITY_OUTPUT_DIR, 'communities_for_analysis.csv'))Rows: 2427 Columns: 12── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): city_name, jetz_species_name, seasonal, presence, origin
dbl (4): city_id, distance_to_northern_edge_km, distance_to_southern_edge_km, relative_abundance_proxy
lgl (3): present_urban_high, present_urban_med, present_urban_low
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.city_introduced_species = communities %>% group_by(city_id) %>% summarise(number_of_species = n()) %>% left_join(
communities %>% group_by(city_id) %>% filter(origin == 'Introduced') %>% summarise(number_of_introduced_species = n())
) %>% replace_na(list(number_of_introduced_species = 0))Joining with `by = join_by(city_id)`community_data_with_introductions = left_join(community_data, city_introduced_species)Joining with `by = join_by(city_id)`community_data_with_introductions$has_introduced_species = community_data_with_introductions$number_of_introduced_species > 0
community_data_with_introductionscommunity_data_with_introductions[,c('mntd_normalised', 'has_introduced_species')]community_data_with_introductions %>% group_by(has_introduced_species) %>% summarise(
total_cities = n(),
mean_mntd_normalised = mean(mntd_normalised, na.rm = T),
median_mntd_normalised = median(mntd_normalised, na.rm = T),
sd_mntd_normalised = sd(mntd_normalised, na.rm = T),
mean_mass_fdiv_normalised = mean(mass_fdiv_normalised, na.rm = T),
median_mass_fdiv_normalised = median(mass_fdiv_normalised, na.rm = T),
sd_mass_fdiv_normalised = sd(mass_fdiv_normalised, na.rm = T),
mean_gape_width_fdiv_normalised = mean(gape_width_fdiv_normalised, na.rm = T),
median_gape_width_fdiv_normalised = median(gape_width_fdiv_normalised, na.rm = T),
sd_gape_width_fdiv_normalised = sd(gape_width_fdiv_normalised, na.rm = T),
mean_handwing_index_fdiv_normalised = mean(handwing_index_fdiv_normalised, na.rm = T),
median_handwing_index_fdiv_normalised = median(handwing_index_fdiv_normalised, na.rm = T),
sd_handwing_index_fdiv_normalised = sd(handwing_index_fdiv_normalised, na.rm = T)
)MNTD
ggplot(community_data_with_introductions, aes(x = has_introduced_species, y = mntd_normalised)) + geom_boxplot()
wilcox.test(mntd_normalised ~ has_introduced_species, community_data_with_introductions, na.action = 'na.omit')
Wilcoxon rank sum test with continuity correction
data: mntd_normalised by has_introduced_species
W = 10512, p-value = 0.02269
alternative hypothesis: true location shift is not equal to 0There is a significant difference between the response of cities with introduced species (0.53±0.27) and those without (0.47±0.19) (p-value = 0.02).
Mass FDiv
ggplot(community_data_with_introductions, aes(x = has_introduced_species, y = mass_fdiv_normalised)) + geom_boxplot()
wilcox.test(mass_fdiv_normalised ~ has_introduced_species, community_data_with_introductions, na.action = 'na.omit')
Wilcoxon rank sum test with continuity correction
data: mass_fdiv_normalised by has_introduced_species
W = 16441, p-value = 0.0000002359
alternative hypothesis: true location shift is not equal to 0There is a significant difference between the response of cities with introduced species (0.57±0.27) and those without (0.73±0.24) (p < 0.0001)
Beak Gape FDiv
ggplot(community_data_with_introductions, aes(x = has_introduced_species, y = gape_width_fdiv_normalised)) + geom_boxplot()
wilcox.test(gape_width_fdiv_normalised ~ has_introduced_species, community_data_with_introductions, na.action = 'na.omit')
Wilcoxon rank sum test with continuity correction
data: gape_width_fdiv_normalised by has_introduced_species
W = 10658, p-value = 0.1538
alternative hypothesis: true location shift is not equal to 0There is NOT a significant difference between the response of cities with introduced species (0.61±0.30) and those without (0.56±0.27)
HWI FDiv
ggplot(community_data_with_introductions, aes(x = has_introduced_species, y = handwing_index_fdiv_normalised)) + geom_boxplot()
wilcox.test(handwing_index_fdiv_normalised ~ has_introduced_species, community_data_with_introductions, na.action = 'na.omit')
Wilcoxon rank sum test with continuity correction
data: handwing_index_fdiv_normalised by has_introduced_species
W = 18856, p-value < 0.00000000000000022
alternative hypothesis: true location shift is not equal to 0There is a significant difference between the response of cities with introduced species (0.49±0.30) and those without (0.79±0.21) (p < 0.0001)
Examine individual metrics
Analysis data frame
geography = read_csv(filename(CITY_DATA_OUTPUT_DIR, 'geography.csv'))Rows: 342 Columns: 26── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
dbl (26): city_id, city_avg_ndvi, city_avg_elevation, city_avg_temp, city_avg_min_monthly_temp, city_avg_max_monthly_temp, city_avg_monthly_temp, city_avg_rai...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.names(geography) [1] "city_id" "city_avg_ndvi" "city_avg_elevation" "city_avg_temp"
[5] "city_avg_min_monthly_temp" "city_avg_max_monthly_temp" "city_avg_monthly_temp" "city_avg_rainfall"
[9] "city_avg_max_monthly_rainfall" "city_avg_min_monthly_rainfall" "city_avg_soil_moisture" "city_max_elev"
[13] "city_min_elev" "city_elev_range" "region_20km_avg_ndvi" "region_20km_avg_elevation"
[17] "region_20km_avg_soil_moisture" "region_20km_max_elev" "region_20km_min_elev" "region_20km_elev_range"
[21] "region_50km_avg_ndvi" "region_50km_avg_elevation" "region_50km_avg_soil_moisture" "region_50km_max_elev"
[25] "region_50km_min_elev" "region_50km_elev_range" analysis_data = community_data_with_realm[,c('city_id', 'mntd_normalised', 'mass_fdiv_normalised', 'gape_width_fdiv_normalised', 'handwing_index_fdiv_normalised', 'core_realm')] %>%
left_join(city_points[,c('city_id', 'latitude', 'longitude')]) %>%
left_join(community_data_with_introductions[,c('city_id', 'has_introduced_species')]) %>%
left_join(geography)Joining with `by = join_by(city_id)`Joining with `by = join_by(city_id)`Joining with `by = join_by(city_id)`analysis_data$abs_latitude = abs(analysis_data$latitude)
analysis_data$core_realm = factor(analysis_data$core_realm, levels = c('Palearctic', 'Nearctic', 'Neotropic', 'Afrotropic', 'Indomalayan', 'Australasia', 'Oceania'))
analysis_data$has_introduced_species = factor(analysis_data$has_introduced_species, level = c('TRUE', 'FALSE'), labels = c('Introduced species', 'No introduced species'))model_data = function(df, dependant_var) {
df[,c(dependant_var, 'core_realm', 'abs_latitude', 'latitude', 'longitude', 'has_introduced_species', 'city_avg_ndvi', 'city_avg_elevation', 'city_avg_temp', 'city_avg_min_monthly_temp', 'city_avg_max_monthly_temp', 'city_avg_monthly_temp', 'city_avg_rainfall', 'city_avg_max_monthly_rainfall', 'city_avg_min_monthly_rainfall', 'city_avg_soil_moisture', 'city_max_elev', 'city_min_elev', 'city_elev_range', 'region_20km_avg_ndvi', 'region_20km_avg_elevation', 'region_20km_avg_soil_moisture', 'region_20km_max_elev', 'region_20km_min_elev', 'region_20km_elev_range', 'region_50km_avg_ndvi', 'region_50km_avg_elevation', 'region_50km_avg_soil_moisture', 'region_50km_max_elev', 'region_50km_min_elev', 'region_50km_elev_range')]
}
model_data(analysis_data, 'mntd_normalised')names(analysis_data) [1] "city_id" "mntd_normalised" "mass_fdiv_normalised" "gape_width_fdiv_normalised"
[5] "handwing_index_fdiv_normalised" "core_realm" "latitude" "longitude"
[9] "geometry" "has_introduced_species" "city_avg_ndvi" "city_avg_elevation"
[13] "city_avg_temp" "city_avg_min_monthly_temp" "city_avg_max_monthly_temp" "city_avg_monthly_temp"
[17] "city_avg_rainfall" "city_avg_max_monthly_rainfall" "city_avg_min_monthly_rainfall" "city_avg_soil_moisture"
[21] "city_max_elev" "city_min_elev" "city_elev_range" "region_20km_avg_ndvi"
[25] "region_20km_avg_elevation" "region_20km_avg_soil_moisture" "region_20km_max_elev" "region_20km_min_elev"
[29] "region_20km_elev_range" "region_50km_avg_ndvi" "region_50km_avg_elevation" "region_50km_avg_soil_moisture"
[33] "region_50km_max_elev" "region_50km_min_elev" "region_50km_elev_range" "abs_latitude" all_explanatories = c(
'abs_latitude', 'latitude', 'longitude',
'has_introduced_species',
'city_avg_ndvi', 'city_avg_elevation', 'city_avg_temp', 'city_avg_min_monthly_temp', 'city_avg_max_monthly_temp',
'city_avg_monthly_temp', 'city_avg_rainfall', 'city_avg_max_monthly_rainfall', 'city_avg_min_monthly_rainfall',
'city_avg_soil_moisture', 'city_max_elev', 'city_min_elev', 'city_elev_range',
'region_20km_avg_ndvi', 'region_20km_avg_elevation', 'region_20km_avg_soil_moisture', 'region_20km_max_elev',
'region_20km_min_elev', 'region_20km_elev_range',
'region_50km_avg_ndvi', 'region_50km_avg_elevation', 'region_50km_avg_soil_moisture', 'region_50km_max_elev',
'region_50km_min_elev', 'region_50km_elev_range',
'core_realmAfrotropic', 'core_realmAustralasia', 'core_realmIndomalayan', 'core_realmNearctic', 'core_realmNeotropic', 'core_realmPalearctic')
type_labels = function(p) {
explanatory_levels = all_explanatories[all_explanatories %in% p$explanatory]
p$explanatory <- factor(p$explanatory, levels = explanatory_levels)
p$type <- 'Realm'
p$type[p$explanatory %in% c('city_avg_ndvi', 'city_avg_elevation', 'city_avg_temp', 'city_avg_min_monthly_temp', 'city_avg_max_monthly_temp',
'city_avg_monthly_temp', 'city_avg_rainfall', 'city_avg_max_monthly_rainfall', 'city_avg_min_monthly_rainfall',
'city_avg_soil_moisture', 'city_max_elev', 'city_min_elev', 'city_elev_range')] <- 'City geography'
p$type[p$explanatory %in% c('region_50km_avg_ndvi', 'region_50km_avg_elevation', 'region_50km_avg_soil_moisture', 'region_50km_max_elev',
'region_50km_min_elev', 'region_50km_elev_range')] <- 'Regional (50km) geography'
p$type[p$explanatory %in% c('region_20km_avg_ndvi', 'region_20km_avg_elevation', 'region_20km_avg_soil_moisture', 'region_20km_max_elev',
'region_20km_min_elev', 'region_20km_elev_range')] <- 'Regional (20km) geography'
p$type[p$explanatory %in% c('abs_latitude', 'latitude', 'longitude')] <- 'Spatial'
p
}explanatory_labels = c(
'has_introduced_species'='Has introduced species',
'city_avg_ndvi'='Average NDVI',
'city_avg_elevation'='Average elevation',
'city_avg_temp'='Average temperature',
'city_avg_min_monthly_temp'='Average minimum monthly temperature',
'city_avg_max_monthly_temp'='Average maximum monthly temperature',
'city_avg_monthly_temp'='Average monthly temperature',
'city_avg_rainfall'='Average rainfall',
'city_avg_max_monthly_rainfall'='Average maximum monthly rainfall',
'city_avg_min_monthly_rainfall'='Average minimum monthly rainfall',
'city_avg_soil_moisture'='Average soil moisture',
'city_max_elev'='Maximum elevation',
'city_min_elev'='Minimum elevation',
'city_elev_range'='Elevation range',
'region_20km_avg_ndvi'='Average NDVI',
'region_20km_avg_elevation'='Average elevation',
'region_20km_avg_soil_moisture'='Average soil moisture',
'region_20km_max_elev'='Maximum elevation',
'region_20km_min_elev'='Minimum elevation',
'region_20km_elev_range'='Elevation range',
'region_50km_avg_ndvi'='Average NDVI',
'region_50km_avg_elevation'='Average elevation',
'region_50km_avg_soil_moisture'='Average soil moisture',
'region_50km_max_elev'='Maximum elevation',
'region_50km_min_elev'='Minimum elevation',
'region_50km_elev_range'='Elevation range',
'abs_latitude' = 'Absolute latitude',
'latitude' = 'Latitude',
'longitude' = 'Longitude',
'core_realmAfrotropic' = 'Afrotropical',
'core_realmAustralasia' = 'Austaliasian',
'core_realmIndomalayan' = 'Indomalayan',
'core_realmNearctic' = 'Nearctic',
'core_realmNeotropic' = 'Neotropical',
'core_realmPalearctic' = 'Palearctic',
'core_realmOceania' = 'Oceanical')Helper plot functions
geom_normalised_histogram = function(name, gg, legend.position = "bottom") {
gg +
geom_histogram(aes(fill = core_realm), binwidth = 0.1, position = "dodge") +
geom_vline(aes(xintercept = 0.5), color = "#000000", size = 0.4) +
geom_vline(aes(xintercept = 0), color = "#000000", size = 0.2, linetype = "dashed") +
geom_vline(aes(xintercept = 1), color = "#000000", size = 0.2, linetype = "dashed") +
ylab("Number of cities") + xlab("Normalised Response") + ylim(c(0, 70)) +
labs(title = name, fill = 'Realm') +
theme_bw() +
theme(legend.position=legend.position)
}geom_map = function(map_sf, title) {
norm_mntd_analysis_geo = ggplot() +
geom_sf(data = world_map, aes(geometry = geometry)) +
map_sf +
normalised_colours_scale +
labs(colour = 'Normalised\nResponse') +
theme_bw() +
theme(legend.position="bottom")
}Helper Dredge functions
# Taken from MuMIN package
# https://rdrr.io/cran/MuMIn/src/R/averaging.R
# https://rdrr.io/cran/MuMIn/src/R/model.avg.R
.coefarr.avg <-
function(cfarr, weight, revised.var, full, alpha) {
weight <- weight / sum(weight)
nCoef <- dim(cfarr)[3L]
if(full) {
nas <- is.na(cfarr[, 1L, ]) & is.na(cfarr[, 2L, ])
cfarr[, 1L, ][nas] <- cfarr[, 2L, ][nas] <- 0
#cfarr[, 1L:2L, ][is.na(cfarr[, 1L:2L, ])] <- 0
if(!all(is.na(cfarr[, 3L, ])))
cfarr[ ,3L, ][is.na(cfarr[ , 3L, ])] <- Inf
}
avgcoef <- array(dim = c(nCoef, 5L),
dimnames = list(dimnames(cfarr)[[3L]], c("Estimate",
"Std. Error", "Adjusted SE", "Lower CI", "Upper CI")))
for(i in seq_len(nCoef))
avgcoef[i, ] <- par.avg(cfarr[, 1L, i], cfarr[, 2L, i], weight,
df = cfarr[, 3L, i], alpha = alpha, revised.var = revised.var)
avgcoef[is.nan(avgcoef)] <- NA
return(avgcoef)
}
.makecoefmat <- function(cf) {
no.ase <- all(is.na(cf[, 3L]))
z <- abs(cf[, 1L] / cf[, if(no.ase) 2L else 3L])
pval <- 2 * pnorm(z, lower.tail = FALSE)
cbind(cf[, if(no.ase) 1L:2L else 1L:3L, drop = FALSE],
`z value` = z, `Pr(>|z|)` = zapsmall(pval))
}
# Generate model selections using lmer, dredge, and model.avg
# `forumla` : a two-sided linear formula object describing both the fixed-effects and random-effects part of the model
# `data` : the data frame containing the variables from the formula
# `aic_delta` : the AIC delta to use for selecting models in model average
model_average <- function(formula, data, aic_delta = 20) {
model <- lm(
formula,
data=data
)
dredge_result <- dredge(model)
summary(model.avg(dredge_result, subset = delta < aic_delta))
}
# Create a summary data frame containing the selected variables from a model
# `model_sum` : The model summary output from `model_average`
model_summary <- function(model_sum) {
.column_name <- function(postfix) {
postfix
}
# just return the estimate and p value
weight <- model_sum$msTable[, 5L]
coefmat.full <- as.data.frame(.makecoefmat(.coefarr.avg(model_sum$coefArray, weight,
attr(model_sum, "revised.var"), TRUE, 0.05)))
coefmat.subset <-
as.data.frame(.makecoefmat(.coefarr.avg(model_sum$coefArray, weight,
attr(model_sum, "revised.var"), FALSE, 0.05)))
coefmat.subset <- coefmat.subset[-c(1), c(1, 2, 5)]
names(coefmat.subset) <- c(.column_name("estimate"), .column_name("error"), .column_name("p"))
coefmat.subset <- tibble::rownames_to_column(coefmat.subset, "explanatory")
coefmat.subset$model = 'subset'
coefmat.full <- coefmat.full[-c(1), c(1, 2, 5)]
names(coefmat.full) <- c(.column_name("estimate"), .column_name("error"), .column_name("p"))
coefmat.full <- tibble::rownames_to_column(coefmat.full, "explanatory")
coefmat.full$model = 'full'
rbind(coefmat.full, coefmat.subset)
}formula_from_vsurp = function(predictors, dependent, vsurp_result) {
as.formula(paste(dependent, paste(names(predictors[,vsurp_result$varselect.interp]), collapse="+"), sep = "~"))
}plot_vsurp_result = function(result_table) {
plot = result_table[result_table$model == 'full',]
plot = type_labels(plot)
ggplot(plot, aes(y=explanatory, x=estimate, colour = type)) +
geom_line() +
geom_point()+
geom_errorbar(aes(xmin=estimate-error, xmax=estimate+error), width=.2,
position=position_dodge(0.05)) +
scale_y_discrete(
limits = rev(levels(plot$explanatory)),
labels = explanatory_labels) +
theme_bw() +
geom_vline(xintercept=0, linetype="dotted") +
guides(colour=guide_legend(title="Predictor type")) + xlab('Increase from 0 (habitat filtering)\nto 1 (competitive exclusion)\n± Standard Error') + ylab('Predictor') +
theme(legend.justification = "top")
}MNTD
norm_mntd_analysis_plot = geom_normalised_histogram(
'MNTD',
ggplot(analysis_data, aes(mntd_normalised))
)Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
Please use `linewidth` instead.norm_mntd_analysis_plot
norm_mntd_analysis_geo_plot = geom_map(geom_sf(data = analysis_data, aes(color = mntd_normalised, geometry = geometry)), 'MNTD')
norm_mntd_analysis_geo_plot
norm_mntd_analysis_data = model_data(analysis_data[!is.na(analysis_data$mntd_normalised),], 'mntd_normalised')
norm_mntd_analysis_predictors = norm_mntd_analysis_data[,-1]
norm_mntd_analysis_interp = VSURF(x = norm_mntd_analysis_predictors, y = norm_mntd_analysis_data$mntd_normalised)Thresholding step
Estimated computational time (on one core): 19.5 sec.
|
| | 0%
|
|======== | 5%
|
|=============== | 10%
|
|======================= | 15%
|
|============================== | 20%
|
|====================================== | 25%
|
|============================================== | 30%
|
|===================================================== | 35%
|
|============================================================= | 40%
|
|==================================================================== | 45%
|
|============================================================================ | 50%
|
|==================================================================================== | 55%
|
|=========================================================================================== | 60%
|
|=================================================================================================== | 65%
|
|========================================================================================================== | 70%
|
|================================================================================================================== | 75%
|
|========================================================================================================================== | 80%
|
|================================================================================================================================= | 85%
|
|========================================================================================================================================= | 90%
|
|================================================================================================================================================ | 95%
|
|========================================================================================================================================================| 100%
Interpretation step (on 30 variables)
Estimated computational time (on one core): between 7.8 sec. and 37.2 sec.
|
| | 0%
|
|===== | 3%
|
|========== | 7%
|
|=============== | 10%
|
|==================== | 13%
|
|========================= | 17%
|
|============================== | 20%
|
|=================================== | 23%
|
|========================================= | 27%
|
|============================================== | 30%
|
|=================================================== | 33%
|
|======================================================== | 37%
|
|============================================================= | 40%
|
|================================================================== | 43%
|
|======================================================================= | 47%
|
|============================================================================ | 50%
|
|================================================================================= | 53%
|
|====================================================================================== | 57%
|
|=========================================================================================== | 60%
|
|================================================================================================ | 63%
|
|===================================================================================================== | 67%
|
|========================================================================================================== | 70%
|
|=============================================================================================================== | 73%
|
|===================================================================================================================== | 77%
|
|========================================================================================================================== | 80%
|
|=============================================================================================================================== | 83%
|
|==================================================================================================================================== | 87%
|
|========================================================================================================================================= | 90%
|
|============================================================================================================================================== | 93%
|
|=================================================================================================================================================== | 97%
|
|========================================================================================================================================================| 100%
Prediction step (on 6 variables)
Maximum estimated computational time (on one core): 2.3 sec.
|
| | 0%
|
|========================= | 17%
|
|=================================================== | 33%
|
|============================================================================ | 50%
|
|===================================================================================================== | 67%
|
|=============================================================================================================================== | 83%
|
|========================================================================================================================================================| 100%names(norm_mntd_analysis_predictors[,norm_mntd_analysis_interp$varselect.interp])[1] "city_avg_monthly_temp" "core_realm" "city_avg_max_monthly_temp" "longitude" "city_avg_temp"
[6] "latitude" norm_mntd_analysis_formula = formula_from_vsurp(norm_mntd_analysis_predictors, "mntd_normalised", norm_mntd_analysis_interp)
norm_mntd_analysis_result <- model_average(norm_mntd_analysis_formula, norm_mntd_analysis_data)Fixed term is "(Intercept)"norm_mntd_analysis_result_table = model_summary(norm_mntd_analysis_result)
norm_mntd_analysis_result_tablenorm_mntd_analysis_pred_plot = plot_vsurp_result(norm_mntd_analysis_result_table)
norm_mntd_analysis_pred_plot
Gape width - FDiv
norm_gape_fdiv_analysis_plot = geom_normalised_histogram(
'Gape FDiv',
ggplot(analysis_data, aes(gape_width_fdiv_normalised))
)
norm_gape_fdiv_analysis_plot
norm_gape_fdiv_analysis_geo_plot = geom_map(geom_sf(data = analysis_data, aes(color = gape_width_fdiv_normalised, geometry = geometry)), 'Gape Width FDiv')
norm_gape_fdiv_analysis_geo_plot
norm_gape_fdiv_analysis_data = model_data(analysis_data[!is.na(analysis_data$gape_width_fdiv_normalised),], 'gape_width_fdiv_normalised')
norm_gape_fdiv_analysis_predictors = norm_gape_fdiv_analysis_data[,-1]
norm_gape_fdiv_analysis_interp = VSURF(x = norm_gape_fdiv_analysis_predictors, y = norm_gape_fdiv_analysis_data$gape_width_fdiv_normalised)Thresholding step
Estimated computational time (on one core): 16.8 sec.
|
| | 0%
|
|======== | 5%
|
|=============== | 10%
|
|======================= | 15%
|
|============================== | 20%
|
|====================================== | 25%
|
|============================================== | 30%
|
|===================================================== | 35%
|
|============================================================= | 40%
|
|==================================================================== | 45%
|
|============================================================================ | 50%
|
|==================================================================================== | 55%
|
|=========================================================================================== | 60%
|
|=================================================================================================== | 65%
|
|========================================================================================================== | 70%
|
|================================================================================================================== | 75%
|
|========================================================================================================================== | 80%
|
|================================================================================================================================= | 85%
|
|========================================================================================================================================= | 90%
|
|================================================================================================================================================ | 95%
|
|========================================================================================================================================================| 100%
Interpretation step (on 30 variables)
Estimated computational time (on one core): between 8.4 sec. and 39 sec.
|
| | 0%
|
|===== | 3%
|
|========== | 7%
|
|=============== | 10%
|
|==================== | 13%
|
|========================= | 17%
|
|============================== | 20%
|
|=================================== | 23%
|
|========================================= | 27%
|
|============================================== | 30%
|
|=================================================== | 33%
|
|======================================================== | 37%
|
|============================================================= | 40%
|
|================================================================== | 43%
|
|======================================================================= | 47%
|
|============================================================================ | 50%
|
|================================================================================= | 53%
|
|====================================================================================== | 57%
|
|=========================================================================================== | 60%
|
|================================================================================================ | 63%
|
|===================================================================================================== | 67%
|
|========================================================================================================== | 70%
|
|=============================================================================================================== | 73%
|
|===================================================================================================================== | 77%
|
|========================================================================================================================== | 80%
|
|=============================================================================================================================== | 83%
|
|==================================================================================================================================== | 87%
|
|========================================================================================================================================= | 90%
|
|============================================================================================================================================== | 93%
|
|=================================================================================================================================================== | 97%
|
|========================================================================================================================================================| 100%
Prediction step (on 9 variables)
Maximum estimated computational time (on one core): 4 sec.
|
| | 0%
|
|================= | 11%
|
|================================== | 22%
|
|=================================================== | 33%
|
|==================================================================== | 44%
|
|==================================================================================== | 56%
|
|===================================================================================================== | 67%
|
|====================================================================================================================== | 78%
|
|======================================================================================================================================= | 89%
|
|========================================================================================================================================================| 100%names(norm_gape_fdiv_analysis_predictors[,norm_gape_fdiv_analysis_interp$varselect.interp])[1] "longitude" "core_realm" "abs_latitude" "city_avg_temp"
[5] "city_avg_max_monthly_temp" "city_avg_min_monthly_temp" "city_avg_max_monthly_rainfall" "latitude"
[9] "city_max_elev" norm_gape_fdiv_analysis_formula = formula_from_vsurp(norm_gape_fdiv_analysis_predictors, "gape_width_fdiv_normalised", norm_gape_fdiv_analysis_interp)
norm_gape_fdiv_analysis_result <- model_average(norm_gape_fdiv_analysis_formula, norm_gape_fdiv_analysis_data)Fixed term is "(Intercept)"norm_gape_fdiv_analysis_result_table = model_summary(norm_gape_fdiv_analysis_result)
norm_gape_fdiv_analysis_result_tablenorm_gape_fdiv_analysis_pred_plot = plot_vsurp_result(norm_gape_fdiv_analysis_result_table)
norm_gape_fdiv_analysis_pred_plot
HWI - FDiv
norm_hwi_fdiv_loco_analysis_plot = geom_normalised_histogram(
'HWI FDiv',
ggplot(analysis_data, aes(handwing_index_fdiv_normalised))
)
norm_hwi_fdiv_loco_analysis_plot
norm_hwi_fdiv_analysis_geo_plot = geom_map(geom_sf(data = analysis_data, aes(color = handwing_index_fdiv_normalised, geometry = geometry)), 'HWI FDiv')
norm_hwi_fdiv_analysis_geo_plot
norm_hwi_fdiv_analysis_data = model_data(analysis_data[!is.na(analysis_data$handwing_index_fdiv_normalised),], 'handwing_index_fdiv_normalised')
norm_hwi_fdiv_analysis_predictors = norm_hwi_fdiv_analysis_data[,-1]
norm_hwi_fdiv_analysis_interp = VSURF(x = norm_hwi_fdiv_analysis_predictors, y = norm_hwi_fdiv_analysis_data$handwing_index_fdiv_normalised)Thresholding step
Estimated computational time (on one core): 15.3 sec.
|
| | 0%
|
|======== | 5%
|
|=============== | 10%
|
|======================= | 15%
|
|============================== | 20%
|
|====================================== | 25%
|
|============================================== | 30%
|
|===================================================== | 35%
|
|============================================================= | 40%
|
|==================================================================== | 45%
|
|============================================================================ | 50%
|
|==================================================================================== | 55%
|
|=========================================================================================== | 60%
|
|=================================================================================================== | 65%
|
|========================================================================================================== | 70%
|
|================================================================================================================== | 75%
|
|========================================================================================================================== | 80%
|
|================================================================================================================================= | 85%
|
|========================================================================================================================================= | 90%
|
|================================================================================================================================================ | 95%
|
|========================================================================================================================================================| 100%
Interpretation step (on 30 variables)
Estimated computational time (on one core): between 6.6 sec. and 36.3 sec.
|
| | 0%
|
|===== | 3%
|
|========== | 7%
|
|=============== | 10%
|
|==================== | 13%
|
|========================= | 17%
|
|============================== | 20%
|
|=================================== | 23%
|
|========================================= | 27%
|
|============================================== | 30%
|
|=================================================== | 33%
|
|======================================================== | 37%
|
|============================================================= | 40%
|
|================================================================== | 43%
|
|======================================================================= | 47%
|
|============================================================================ | 50%
|
|================================================================================= | 53%
|
|====================================================================================== | 57%
|
|=========================================================================================== | 60%
|
|================================================================================================ | 63%
|
|===================================================================================================== | 67%
|
|========================================================================================================== | 70%
|
|=============================================================================================================== | 73%
|
|===================================================================================================================== | 77%
|
|========================================================================================================================== | 80%
|
|=============================================================================================================================== | 83%
|
|==================================================================================================================================== | 87%
|
|========================================================================================================================================= | 90%
|
|============================================================================================================================================== | 93%
|
|=================================================================================================================================================== | 97%
|
|========================================================================================================================================================| 100%
Prediction step (on 7 variables)
Maximum estimated computational time (on one core): 2.3 sec.
|
| | 0%
|
|====================== | 14%
|
|=========================================== | 29%
|
|================================================================= | 43%
|
|======================================================================================= | 57%
|
|============================================================================================================= | 71%
|
|================================================================================================================================== | 86%
|
|========================================================================================================================================================| 100%names(norm_hwi_fdiv_analysis_predictors[,norm_hwi_fdiv_analysis_interp$varselect.interp])[1] "latitude" "core_realm" "city_avg_max_monthly_temp" "city_avg_temp"
[5] "longitude" "city_avg_max_monthly_rainfall" "region_50km_min_elev" norm_hwi_fdiv_analysis_formula = formula_from_vsurp(norm_hwi_fdiv_analysis_predictors, "handwing_index_fdiv_normalised", norm_hwi_fdiv_analysis_interp)
norm_hwi_fdiv_analysis_result <- model_average(norm_hwi_fdiv_analysis_formula, norm_hwi_fdiv_analysis_data)Fixed term is "(Intercept)"norm_hwi_fdiv_analysis_result_table = model_summary(norm_hwi_fdiv_analysis_result)
norm_hwi_fdiv_analysis_result_tablenorm_hwi_fdiv_analysis_pred_plot = plot_vsurp_result(norm_hwi_fdiv_analysis_result_table)
norm_hwi_fdiv_analysis_pred_plot
Mass - FDiv
norm_mass_fdiv_loco_analysis_plot = geom_normalised_histogram(
'Mass FDiv',
ggplot(analysis_data, aes(mass_fdiv_normalised))
)
norm_mass_fdiv_loco_analysis_plot
norm_mass_fdiv_analysis_geo_plot = geom_map(geom_sf(data = analysis_data, aes(color = mass_fdiv_normalised, geometry = geometry)), 'Mass FDiv')
norm_mass_fdiv_analysis_geo_plot
norm_mass_fdiv_analysis_data = model_data(analysis_data[!is.na(analysis_data$mass_fdiv_normalised),], 'mass_fdiv_normalised')
norm_mass_fdiv_analysis_predictors = norm_mass_fdiv_analysis_data[,-1]
norm_mass_fdiv_analysis_interp = VSURF(x = norm_mass_fdiv_analysis_predictors, y = norm_mass_fdiv_analysis_data$mass_fdiv_normalised)Thresholding step
Estimated computational time (on one core): 16.6 sec.
|
| | 0%
|
|======== | 5%
|
|=============== | 10%
|
|======================= | 15%
|
|============================== | 20%
|
|====================================== | 25%
|
|============================================== | 30%
|
|===================================================== | 35%
|
|============================================================= | 40%
|
|==================================================================== | 45%
|
|============================================================================ | 50%
|
|==================================================================================== | 55%
|
|=========================================================================================== | 60%
|
|=================================================================================================== | 65%
|
|========================================================================================================== | 70%
|
|================================================================================================================== | 75%
|
|========================================================================================================================== | 80%
|
|================================================================================================================================= | 85%
|
|========================================================================================================================================= | 90%
|
|================================================================================================================================================ | 95%
|
|========================================================================================================================================================| 100%
Interpretation step (on 30 variables)
Estimated computational time (on one core): between 1.8 sec. and 34.5 sec.
|
| | 0%
|
|===== | 3%
|
|========== | 7%
|
|=============== | 10%
|
|==================== | 13%
|
|========================= | 17%
|
|============================== | 20%
|
|=================================== | 23%
|
|========================================= | 27%
|
|============================================== | 30%
|
|=================================================== | 33%
|
|======================================================== | 37%
|
|============================================================= | 40%
|
|================================================================== | 43%
|
|======================================================================= | 47%
|
|============================================================================ | 50%
|
|================================================================================= | 53%
|
|====================================================================================== | 57%
|
|=========================================================================================== | 60%
|
|================================================================================================ | 63%
|
|===================================================================================================== | 67%
|
|========================================================================================================== | 70%
|
|=============================================================================================================== | 73%
|
|===================================================================================================================== | 77%
|
|========================================================================================================================== | 80%
|
|=============================================================================================================================== | 83%
|
|==================================================================================================================================== | 87%
|
|========================================================================================================================================= | 90%
|
|============================================================================================================================================== | 93%
|
|=================================================================================================================================================== | 97%
|
|========================================================================================================================================================| 100%
Prediction step (on 5 variables)
Maximum estimated computational time (on one core): 0.9 sec.
|
| | 0%
|
|============================== | 20%
|
|============================================================= | 40%
|
|=========================================================================================== | 60%
|
|========================================================================================================================== | 80%
|
|========================================================================================================================================================| 100%names(norm_mass_fdiv_analysis_predictors[,norm_mass_fdiv_analysis_interp$varselect.interp])[1] "core_realm" "abs_latitude" "city_avg_temp" "latitude"
[5] "city_avg_max_monthly_rainfall"norm_mass_fdiv_analysis_formula = formula_from_vsurp(norm_mass_fdiv_analysis_predictors, "mass_fdiv_normalised", norm_mass_fdiv_analysis_interp)
norm_mass_fdiv_analysis_result <- model_average(norm_mass_fdiv_analysis_formula, norm_mass_fdiv_analysis_data)Fixed term is "(Intercept)"norm_mass_fdiv_analysis_result_table = model_summary(norm_mass_fdiv_analysis_result)
norm_mass_fdiv_analysis_result_tablenorm_mass_fdiv_analysis_pred_plot = plot_vsurp_result(norm_mass_fdiv_analysis_result_table)
norm_mass_fdiv_analysis_pred_plot
Create plot of differences in process response
pred_legend <- get_legend(
# create some space to the left of the legend
norm_hwi_fdiv_analysis_pred_plot + theme(legend.box.margin = margin(0, 0, 0, 12))
)Warning: Removed 1 row containing missing values or values outside the scale range (`geom_line()`).`geom_line()`: Each group consists of only one observation.
ℹ Do you need to adjust the group aesthetic?Warning: Removed 1 row containing missing values or values outside the scale range (`geom_point()`).Warning: Multiple components found; returning the first one. To return all, use `return_all = TRUE`.geo_legend <- get_legend(
# create some space to the left of the legend
norm_mass_fdiv_analysis_geo_plot + theme(legend.box.margin = margin(0, 0, 0, 12))
)Warning: Multiple components found; returning the first one. To return all, use `return_all = TRUE`.plot_grid(
plot_grid(
norm_mntd_analysis_geo_plot + theme(legend.position="none"),
norm_mntd_analysis_pred_plot + theme(legend.position="none"),
nrow = 1
) + draw_label("MNTD", size = 16, angle = 90, x = 0.01, y = 0.5),
plot_grid(
norm_gape_fdiv_analysis_geo_plot + theme(legend.position="none"),
norm_gape_fdiv_analysis_pred_plot + theme(legend.position="none"),
nrow = 1
) + draw_label("Gape", size = 16, angle = 90, x = 0.01, y = 0.5),
plot_grid(
norm_hwi_fdiv_analysis_geo_plot + theme(legend.position="none"),
norm_hwi_fdiv_analysis_pred_plot + theme(legend.position="none"),
nrow = 1
) + draw_label("HWI", size = 16, angle = 90, x = 0.01, y = 0.5),
plot_grid(
norm_mass_fdiv_analysis_geo_plot + theme(legend.position="none"),
norm_mass_fdiv_analysis_pred_plot + theme(legend.position="none"),
nrow = 1
) + draw_label("Mass", size = 16, angle = 90, x = 0.01, y = 0.5),
plot_grid(
geo_legend,
pred_legend,
nrow = 1
),
nrow = 5
)`geom_line()`: Each group consists of only one observation.
ℹ Do you need to adjust the group aesthetic?`geom_line()`: Each group consists of only one observation.
ℹ Do you need to adjust the group aesthetic?Warning: Removed 1 row containing missing values or values outside the scale range (`geom_line()`).`geom_line()`: Each group consists of only one observation.
ℹ Do you need to adjust the group aesthetic?Warning: Removed 1 row containing missing values or values outside the scale range (`geom_point()`).Warning: Removed 1 row containing missing values or values outside the scale range (`geom_line()`).`geom_line()`: Each group consists of only one observation.
ℹ Do you need to adjust the group aesthetic?Warning: Removed 1 row containing missing values or values outside the scale range (`geom_point()`).ggsave(filename(FIGURES_OUTPUT_DIR, 'process_response.jpg'), width = 2500, height=5000, units = 'px')
Compare metrics against each other
ggplot(analysis_data, aes(x = gape_width_fdiv_normalised, y = mntd_normalised, colour = core_realm)) +
geom_point() +
ylab("MNTD") +
xlab("Gape Width FDiv") +
theme_bw() + labs(color = "Realm")
ggplot(analysis_data, aes(x = handwing_index_fdiv_normalised, y = mntd_normalised, colour = core_realm)) +
geom_point() +
ylab("MNTD") +
xlab("HWI FDiv") +
theme_bw() + labs(color = "Realm")
ggplot(analysis_data, aes(x = handwing_index_fdiv_normalised, y = gape_width_fdiv_normalised, colour = core_realm)) +
geom_point() +
ylab("Gape Width FDiv") +
xlab("HWI FDiv") +
theme_bw() + labs(color = "Realm")
mntd_fdiv_analysis = analysis_data %>%
dplyr::select(city_id, mntd_normalised, handwing_index_fdiv_normalised, gape_width_fdiv_normalised) %>%
left_join(community_summary) %>%
mutate(urban_pool_perc = urban_pool_size * 100 / regional_pool_size)Joining with `by = join_by(city_id)`mntd_fdiv_analysisggpairs(mntd_fdiv_analysis %>% dplyr::select(mntd_normalised, handwing_index_fdiv_normalised, gape_width_fdiv_normalised, regional_pool_size, urban_pool_size, urban_pool_perc))
ggsave(filename(FIGURES_OUTPUT_DIR, 'appendix_normalised_correlation.jpg'))Saving 7.29 x 4.51 in image
LS0tCnRpdGxlOiAiTWV0cmljcyBmb3IgYXNzZXNzaW5nIGNvbW11bml0eSBhc3NlbWJseSBwcm9jZXNzZXMiCm91dHB1dDogaHRtbF9ub3RlYm9vawpiaWJsaW9ncmFwaHk6IC4uL3JlZi5iaWIgIAotLS0KCmBgYHtyfQpzb3VyY2UoJy4uL2Vudi5SJykKYGBgCgpgYGB7cn0KY29tbXVuaXR5X2RhdGEgPSByZWFkX2NzdihmaWxlbmFtZShDT01NVU5JVFlfT1VUUFVUX0RJUiwgJ2NvbW11bml0eV9hc3NlbWJseV9tZXRyaWNzX3VzaW5nX3JlbGF0aXZlX2FidW5kYW5jZS5jc3YnKSkKaGVhZChjb21tdW5pdHlfZGF0YSkKY29sbmFtZXMoY29tbXVuaXR5X2RhdGEpCmBgYAoKSm9pbiBvbiByZWFsbXMKYGBge3J9CmNpdHlfdG9fcmVhbG0gPSByZWFkX2NzdihmaWxlbmFtZShDSVRZX0RBVEFfT1VUUFVUX0RJUiwgJ3JlYWxtcy5jc3YnKSkKY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSA9IGxlZnRfam9pbihjb21tdW5pdHlfZGF0YSwgY2l0eV90b19yZWFsbSkKYGBgCgpDaXRpZXMgYXMgcG9pbnRzCmBgYHtyfQpjaXR5X3BvaW50cyA9IHN0X2NlbnRyb2lkKHJlYWRfc2YoZmlsZW5hbWUoQ0lUWV9EQVRBX09VVFBVVF9ESVIsICdjaXR5X3NlbGVjdGlvbi5zaHAnKSkpICU+JSBsZWZ0X2pvaW4oY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSkKY2l0eV9wb2ludHNfY29vcmRzID0gc3RfY29vcmRpbmF0ZXMoY2l0eV9wb2ludHMpCmNpdHlfcG9pbnRzJGxhdGl0dWRlID0gY2l0eV9wb2ludHNfY29vcmRzWywxXQpjaXR5X3BvaW50cyRsb25naXR1ZGUgPSBjaXR5X3BvaW50c19jb29yZHNbLDJdCmBgYAogIApgYGB7cn0Kd29ybGRfbWFwID0gcmVhZF9jb3VudHJ5X2JvdW5kYXJpZXMoKQpgYGAKCkxvYWQgY29tbXVuaXR5IGRhdGEsIGFuZCBjcmVhdGUgbG9uZyBmb3JtYXQgdmVyc2lvbgpgYGB7cn0KY29tbXVuaXRpZXMgPSByZWFkX2NzdihmaWxlbmFtZShDT01NVU5JVFlfT1VUUFVUX0RJUiwgJ2NvbW11bml0aWVzX2Zvcl9hbmFseXNpcy5jc3YnKSkKY29tbXVuaXRpZXMKYGBgCgpgYGB7cn0KY29tbXVuaXR5X3N1bW1hcnkgPSBjb21tdW5pdGllcyAlPiUgZ3JvdXBfYnkoY2l0eV9pZCkgJT4lIHN1bW1hcmlzZShyZWdpb25hbF9wb29sX3NpemUgPSBuKCksIHVyYmFuX3Bvb2xfc2l6ZSA9IHN1bShyZWxhdGl2ZV9hYnVuZGFuY2VfcHJveHkgPiAwKSkKY29tbXVuaXR5X3N1bW1hcnkKYGBgCgpMb2FkIHRyYWl0IGRhdGEKYGBge3J9CnRyYWl0cyA9IHJlYWRfY3N2KGZpbGVuYW1lKFRBWE9OT01ZX09VVFBVVF9ESVIsICd0cmFpdHNfamV0ei5jc3YnKSkKaGVhZCh0cmFpdHMpCmBgYAoKTG9hZCByZWFsbSBnZW8KYGBge3J9CnJlc29sdmUgPSByZWFkX3Jlc29sdmUoKQpoZWFkKHJlc29sdmUpCmBgYAoKIyBTdW1tYXJ5IG1ldHJpY3MgYnkgUmVhbG0KYGBge3J9CnRlc3RfcmVxdWlyZWRfdmFsdWVzID0gZnVuY3Rpb24obmFtZSwgZGYpIHsKICBjYXQocGFzdGUoCiAgICB0ZXN0X3ZhbHVlX3dpbGNveChwYXN0ZShuYW1lLCAnTU5URCcpLCBkZiRtbnRkX25vcm1hbGlzZWQpLAogICAgdGVzdF92YWx1ZV93aWxjb3gocGFzdGUobmFtZSwgJ0JlYWsgR2FwZSBGRGl2JyksIGRmJGdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkKSwKICAgIHRlc3RfdmFsdWVfd2lsY294KHBhc3RlKG5hbWUsICdIV0kgRkRpdicpLCBkZiRoYW5kd2luZ19pbmRleF9mZGl2X25vcm1hbGlzZWQpLAogICAgdGVzdF92YWx1ZV93aWxjb3gocGFzdGUobmFtZSwgJ01hc3MgRkRpdicpLCBkZiRtYXNzX2ZkaXZfbm9ybWFsaXNlZCksCiAgICBucm93KGRmKSwKICAgIHNlcCA9ICJcbiIpKQp9CmBgYAoKYGBge3J9CnRlc3RfcmVxdWlyZWRfdmFsdWVzKCdHbG9iYWwnLCBjb21tdW5pdHlfZGF0YV93aXRoX3JlYWxtKQpgYGAKCmBgYHtyfQp1bmlxdWUoY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSRjb3JlX3JlYWxtKQpgYGAKCmBgYHtyfQp0ZXN0X3JlcXVpcmVkX3ZhbHVlcygnTmVhcmN0aWMnLCBjb21tdW5pdHlfZGF0YV93aXRoX3JlYWxtW2NvbW11bml0eV9kYXRhX3dpdGhfcmVhbG0kY29yZV9yZWFsbSA9PSAnTmVhcmN0aWMnLF0pCmBgYAoKYGBge3J9CnRlc3RfcmVxdWlyZWRfdmFsdWVzKCdOZW90cm9waWMnLCBjb21tdW5pdHlfZGF0YV93aXRoX3JlYWxtW2NvbW11bml0eV9kYXRhX3dpdGhfcmVhbG0kY29yZV9yZWFsbSA9PSAnTmVvdHJvcGljJyxdKQpgYGAKCmBgYHtyfQp0ZXN0X3JlcXVpcmVkX3ZhbHVlcygnUGFsZWFyY3RpYycsIGNvbW11bml0eV9kYXRhX3dpdGhfcmVhbG1bY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSRjb3JlX3JlYWxtID09ICdQYWxlYXJjdGljJyxdKQpgYGAKCmBgYHtyfQp0ZXN0X3JlcXVpcmVkX3ZhbHVlcygnQWZyb3Ryb3BpYycsIGNvbW11bml0eV9kYXRhX3dpdGhfcmVhbG1bY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSRjb3JlX3JlYWxtID09ICdBZnJvdHJvcGljJyxdKQpgYGAKCmBgYHtyfQp0ZXN0X3JlcXVpcmVkX3ZhbHVlcygnSW5kb21hbGF5YW4nLCBjb21tdW5pdHlfZGF0YV93aXRoX3JlYWxtW2NvbW11bml0eV9kYXRhX3dpdGhfcmVhbG0kY29yZV9yZWFsbSA9PSAnSW5kb21hbGF5YW4nLF0pCmBgYAoKYGBge3J9CnRlc3RfcmVxdWlyZWRfdmFsdWVzKCdBdXN0cmFsYXNpYScsIGNvbW11bml0eV9kYXRhX3dpdGhfcmVhbG1bY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbSRjb3JlX3JlYWxtID09ICdBdXN0cmFsYXNpYScsXSkKYGBgCgojIFdoYXQgZmFtaWxpZXMgZXhpc3QgaW4gd2hpY2ggcmVhbG1zPwpgYGB7cn0KY29tbXVuaXRpZXMgJT4lIAogIGxlZnRfam9pbihjaXR5X3RvX3JlYWxtKSAlPiUgCiAgbXV0YXRlKGZhbWlseSA9IGdzdWIoICJfLiokIiwgIiIsIGpldHpfc3BlY2llc19uYW1lKSkgJT4lCiAgZHBseXI6OnNlbGVjdChmYW1pbHksIGNvcmVfcmVhbG0pICU+JQogIGRpc3RpbmN0KCkgJT4lCiAgYXJyYW5nZShjb3JlX3JlYWxtKQpgYGAKCiMgU3VtbWFyeSBtZXRyaWNzIGJ5IGludHJvZHVjZWQgc3BlY2llcwpgYGB7cn0KY29tbXVuaXRpZXMgPSByZWFkX2NzdihmaWxlbmFtZShDT01NVU5JVFlfT1VUUFVUX0RJUiwgJ2NvbW11bml0aWVzX2Zvcl9hbmFseXNpcy5jc3YnKSkKY2l0eV9pbnRyb2R1Y2VkX3NwZWNpZXMgPSBjb21tdW5pdGllcyAlPiUgZ3JvdXBfYnkoY2l0eV9pZCkgJT4lIHN1bW1hcmlzZShudW1iZXJfb2Zfc3BlY2llcyA9IG4oKSkgJT4lIGxlZnRfam9pbigKICBjb21tdW5pdGllcyAlPiUgZ3JvdXBfYnkoY2l0eV9pZCkgJT4lIGZpbHRlcihvcmlnaW4gPT0gJ0ludHJvZHVjZWQnKSAlPiUgc3VtbWFyaXNlKG51bWJlcl9vZl9pbnRyb2R1Y2VkX3NwZWNpZXMgPSBuKCkpCikgJT4lIHJlcGxhY2VfbmEobGlzdChudW1iZXJfb2ZfaW50cm9kdWNlZF9zcGVjaWVzID0gMCkpCgpjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMgPSBsZWZ0X2pvaW4oY29tbXVuaXR5X2RhdGEsIGNpdHlfaW50cm9kdWNlZF9zcGVjaWVzKQpjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMkaGFzX2ludHJvZHVjZWRfc3BlY2llcyA9IGNvbW11bml0eV9kYXRhX3dpdGhfaW50cm9kdWN0aW9ucyRudW1iZXJfb2ZfaW50cm9kdWNlZF9zcGVjaWVzID4gMApjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMKYGBgCgpgYGB7cn0KY29tbXVuaXR5X2RhdGFfd2l0aF9pbnRyb2R1Y3Rpb25zWyxjKCdtbnRkX25vcm1hbGlzZWQnLCAnaGFzX2ludHJvZHVjZWRfc3BlY2llcycpXQpgYGAKCmBgYHtyfQpjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMgJT4lIGdyb3VwX2J5KGhhc19pbnRyb2R1Y2VkX3NwZWNpZXMpICU+JSBzdW1tYXJpc2UoCiAgdG90YWxfY2l0aWVzID0gbigpLCAKICAKICBtZWFuX21udGRfbm9ybWFsaXNlZCA9IG1lYW4obW50ZF9ub3JtYWxpc2VkLCBuYS5ybSA9IFQpLAogIG1lZGlhbl9tbnRkX25vcm1hbGlzZWQgPSBtZWRpYW4obW50ZF9ub3JtYWxpc2VkLCBuYS5ybSA9IFQpLAogIHNkX21udGRfbm9ybWFsaXNlZCA9IHNkKG1udGRfbm9ybWFsaXNlZCwgbmEucm0gPSBUKSwKICAKICBtZWFuX21hc3NfZmRpdl9ub3JtYWxpc2VkID0gbWVhbihtYXNzX2ZkaXZfbm9ybWFsaXNlZCwgbmEucm0gPSBUKSwKICBtZWRpYW5fbWFzc19mZGl2X25vcm1hbGlzZWQgPSBtZWRpYW4obWFzc19mZGl2X25vcm1hbGlzZWQsIG5hLnJtID0gVCksCiAgc2RfbWFzc19mZGl2X25vcm1hbGlzZWQgPSBzZChtYXNzX2ZkaXZfbm9ybWFsaXNlZCwgbmEucm0gPSBUKSwKICAKICBtZWFuX2dhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkID0gbWVhbihnYXBlX3dpZHRoX2ZkaXZfbm9ybWFsaXNlZCwgbmEucm0gPSBUKSwKICBtZWRpYW5fZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQgPSBtZWRpYW4oZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQsIG5hLnJtID0gVCksCiAgc2RfZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQgPSBzZChnYXBlX3dpZHRoX2ZkaXZfbm9ybWFsaXNlZCwgbmEucm0gPSBUKSwKICAKICBtZWFuX2hhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCA9IG1lYW4oaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkLCBuYS5ybSA9IFQpLAogIG1lZGlhbl9oYW5kd2luZ19pbmRleF9mZGl2X25vcm1hbGlzZWQgPSBtZWRpYW4oaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkLCBuYS5ybSA9IFQpLAogIHNkX2hhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCA9IHNkKGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCwgbmEucm0gPSBUKQopCmBgYAoKIyMgTU5URApgYGB7cn0KZ2dwbG90KGNvbW11bml0eV9kYXRhX3dpdGhfaW50cm9kdWN0aW9ucywgYWVzKHggPSBoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCB5ID0gbW50ZF9ub3JtYWxpc2VkKSkgKyBnZW9tX2JveHBsb3QoKQpgYGAKCmBgYHtyfQp3aWxjb3gudGVzdChtbnRkX25vcm1hbGlzZWQgfiBoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCBjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMsIG5hLmFjdGlvbiA9ICduYS5vbWl0JykKYGBgCgpUaGVyZSBpcyBhIHNpZ25pZmljYW50IGRpZmZlcmVuY2UgYmV0d2VlbiB0aGUgcmVzcG9uc2Ugb2YgY2l0aWVzIHdpdGggaW50cm9kdWNlZCBzcGVjaWVzICgwLjUzwrEwLjI3KSBhbmQgdGhvc2Ugd2l0aG91dCAoMC40N8KxMC4xOSkgKHAtdmFsdWUgPSAwLjAyKS4KCgojIyBNYXNzIEZEaXYKYGBge3J9CmdncGxvdChjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMsIGFlcyh4ID0gaGFzX2ludHJvZHVjZWRfc3BlY2llcywgeSA9IG1hc3NfZmRpdl9ub3JtYWxpc2VkKSkgKyBnZW9tX2JveHBsb3QoKQpgYGAKCmBgYHtyfQp3aWxjb3gudGVzdChtYXNzX2ZkaXZfbm9ybWFsaXNlZCB+IGhhc19pbnRyb2R1Y2VkX3NwZWNpZXMsIGNvbW11bml0eV9kYXRhX3dpdGhfaW50cm9kdWN0aW9ucywgbmEuYWN0aW9uID0gJ25hLm9taXQnKQpgYGAKVGhlcmUgaXMgYSBzaWduaWZpY2FudCBkaWZmZXJlbmNlIGJldHdlZW4gdGhlIHJlc3BvbnNlIG9mIGNpdGllcyB3aXRoIGludHJvZHVjZWQgc3BlY2llcyAoMC41N8KxMC4yNykgYW5kIHRob3NlIHdpdGhvdXQgKDAuNzPCsTAuMjQpIChwIDwgMC4wMDAxKQoKCiMjIEJlYWsgR2FwZSBGRGl2CmBgYHtyfQpnZ3Bsb3QoY29tbXVuaXR5X2RhdGFfd2l0aF9pbnRyb2R1Y3Rpb25zLCBhZXMoeCA9IGhhc19pbnRyb2R1Y2VkX3NwZWNpZXMsIHkgPSBnYXBlX3dpZHRoX2ZkaXZfbm9ybWFsaXNlZCkpICsgZ2VvbV9ib3hwbG90KCkKYGBgCgpgYGB7cn0Kd2lsY294LnRlc3QoZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQgfiBoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCBjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMsIG5hLmFjdGlvbiA9ICduYS5vbWl0JykKYGBgClRoZXJlIGlzIE5PVCBhIHNpZ25pZmljYW50IGRpZmZlcmVuY2UgYmV0d2VlbiB0aGUgcmVzcG9uc2Ugb2YgY2l0aWVzIHdpdGggaW50cm9kdWNlZCBzcGVjaWVzICgwLjYxwrEwLjMwKSBhbmQgdGhvc2Ugd2l0aG91dCAoMC41NsKxMC4yNykKCgojIyBIV0kgRkRpdgpgYGB7cn0KZ2dwbG90KGNvbW11bml0eV9kYXRhX3dpdGhfaW50cm9kdWN0aW9ucywgYWVzKHggPSBoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCB5ID0gaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkKSkgKyBnZW9tX2JveHBsb3QoKQpgYGAKCmBgYHtyfQp3aWxjb3gudGVzdChoYW5kd2luZ19pbmRleF9mZGl2X25vcm1hbGlzZWQgfiBoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCBjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnMsIG5hLmFjdGlvbiA9ICduYS5vbWl0JykKYGBgClRoZXJlIGlzIGEgc2lnbmlmaWNhbnQgZGlmZmVyZW5jZSBiZXR3ZWVuIHRoZSByZXNwb25zZSBvZiBjaXRpZXMgd2l0aCBpbnRyb2R1Y2VkIHNwZWNpZXMgKDAuNDnCsTAuMzApIGFuZCB0aG9zZSB3aXRob3V0ICgwLjc5wrEwLjIxKSAocCA8IDAuMDAwMSkKCgojIEV4YW1pbmUgaW5kaXZpZHVhbCBtZXRyaWNzCgojIyBBbmFseXNpcyBkYXRhIGZyYW1lCmBgYHtyfQpnZW9ncmFwaHkgPSByZWFkX2NzdihmaWxlbmFtZShDSVRZX0RBVEFfT1VUUFVUX0RJUiwgJ2dlb2dyYXBoeS5jc3YnKSkKbmFtZXMoZ2VvZ3JhcGh5KQpgYGAKCmBgYHtyfQphbmFseXNpc19kYXRhID0gY29tbXVuaXR5X2RhdGFfd2l0aF9yZWFsbVssYygnY2l0eV9pZCcsICdtbnRkX25vcm1hbGlzZWQnLCAnbWFzc19mZGl2X25vcm1hbGlzZWQnLCAnZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQnLCAnaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkJywgJ2NvcmVfcmVhbG0nKV0gJT4lIAogIGxlZnRfam9pbihjaXR5X3BvaW50c1ssYygnY2l0eV9pZCcsICdsYXRpdHVkZScsICdsb25naXR1ZGUnKV0pICU+JQogIGxlZnRfam9pbihjb21tdW5pdHlfZGF0YV93aXRoX2ludHJvZHVjdGlvbnNbLGMoJ2NpdHlfaWQnLCAnaGFzX2ludHJvZHVjZWRfc3BlY2llcycpXSkgJT4lCiAgbGVmdF9qb2luKGdlb2dyYXBoeSkKCmFuYWx5c2lzX2RhdGEkYWJzX2xhdGl0dWRlID0gYWJzKGFuYWx5c2lzX2RhdGEkbGF0aXR1ZGUpCmFuYWx5c2lzX2RhdGEkY29yZV9yZWFsbSA9IGZhY3RvcihhbmFseXNpc19kYXRhJGNvcmVfcmVhbG0sIGxldmVscyA9IGMoJ1BhbGVhcmN0aWMnLCAnTmVhcmN0aWMnLCAnTmVvdHJvcGljJywgJ0Fmcm90cm9waWMnLCAnSW5kb21hbGF5YW4nLCAnQXVzdHJhbGFzaWEnLCAnT2NlYW5pYScpKQphbmFseXNpc19kYXRhJGhhc19pbnRyb2R1Y2VkX3NwZWNpZXMgPSBmYWN0b3IoYW5hbHlzaXNfZGF0YSRoYXNfaW50cm9kdWNlZF9zcGVjaWVzLCBsZXZlbCA9IGMoJ1RSVUUnLCAnRkFMU0UnKSwgbGFiZWxzID0gYygnSW50cm9kdWNlZCBzcGVjaWVzJywgJ05vIGludHJvZHVjZWQgc3BlY2llcycpKQpgYGAKCmBgYHtyfQptb2RlbF9kYXRhID0gZnVuY3Rpb24oZGYsIGRlcGVuZGFudF92YXIpIHsKICBkZlssYyhkZXBlbmRhbnRfdmFyLCAnY29yZV9yZWFsbScsICdhYnNfbGF0aXR1ZGUnLCAnbGF0aXR1ZGUnLCAnbG9uZ2l0dWRlJywgJ2hhc19pbnRyb2R1Y2VkX3NwZWNpZXMnLCAnY2l0eV9hdmdfbmR2aScsICdjaXR5X2F2Z19lbGV2YXRpb24nLCAnY2l0eV9hdmdfdGVtcCcsICdjaXR5X2F2Z19taW5fbW9udGhseV90ZW1wJywgJ2NpdHlfYXZnX21heF9tb250aGx5X3RlbXAnLCAnY2l0eV9hdmdfbW9udGhseV90ZW1wJywgJ2NpdHlfYXZnX3JhaW5mYWxsJywgJ2NpdHlfYXZnX21heF9tb250aGx5X3JhaW5mYWxsJywgJ2NpdHlfYXZnX21pbl9tb250aGx5X3JhaW5mYWxsJywgJ2NpdHlfYXZnX3NvaWxfbW9pc3R1cmUnLCAnY2l0eV9tYXhfZWxldicsICdjaXR5X21pbl9lbGV2JywgJ2NpdHlfZWxldl9yYW5nZScsICdyZWdpb25fMjBrbV9hdmdfbmR2aScsICdyZWdpb25fMjBrbV9hdmdfZWxldmF0aW9uJywgJ3JlZ2lvbl8yMGttX2F2Z19zb2lsX21vaXN0dXJlJywgJ3JlZ2lvbl8yMGttX21heF9lbGV2JywgJ3JlZ2lvbl8yMGttX21pbl9lbGV2JywgJ3JlZ2lvbl8yMGttX2VsZXZfcmFuZ2UnLCAncmVnaW9uXzUwa21fYXZnX25kdmknLCAncmVnaW9uXzUwa21fYXZnX2VsZXZhdGlvbicsICdyZWdpb25fNTBrbV9hdmdfc29pbF9tb2lzdHVyZScsICdyZWdpb25fNTBrbV9tYXhfZWxldicsICdyZWdpb25fNTBrbV9taW5fZWxldicsICdyZWdpb25fNTBrbV9lbGV2X3JhbmdlJyldCn0KbW9kZWxfZGF0YShhbmFseXNpc19kYXRhLCAnbW50ZF9ub3JtYWxpc2VkJykKYGBgCgpgYGB7cn0KbmFtZXMoYW5hbHlzaXNfZGF0YSkKYGBgCgpgYGB7cn0KYWxsX2V4cGxhbmF0b3JpZXMgPSBjKAogICAgJ2Fic19sYXRpdHVkZScsICdsYXRpdHVkZScsICdsb25naXR1ZGUnLCAKICAgICdoYXNfaW50cm9kdWNlZF9zcGVjaWVzJywKICAgICdjaXR5X2F2Z19uZHZpJywgJ2NpdHlfYXZnX2VsZXZhdGlvbicsICdjaXR5X2F2Z190ZW1wJywgJ2NpdHlfYXZnX21pbl9tb250aGx5X3RlbXAnLCAnY2l0eV9hdmdfbWF4X21vbnRobHlfdGVtcCcsIAogICAgJ2NpdHlfYXZnX21vbnRobHlfdGVtcCcsICdjaXR5X2F2Z19yYWluZmFsbCcsICdjaXR5X2F2Z19tYXhfbW9udGhseV9yYWluZmFsbCcsICdjaXR5X2F2Z19taW5fbW9udGhseV9yYWluZmFsbCcsIAogICAgJ2NpdHlfYXZnX3NvaWxfbW9pc3R1cmUnLCAnY2l0eV9tYXhfZWxldicsICdjaXR5X21pbl9lbGV2JywgJ2NpdHlfZWxldl9yYW5nZScsCiAgICAncmVnaW9uXzIwa21fYXZnX25kdmknLCAncmVnaW9uXzIwa21fYXZnX2VsZXZhdGlvbicsICdyZWdpb25fMjBrbV9hdmdfc29pbF9tb2lzdHVyZScsICdyZWdpb25fMjBrbV9tYXhfZWxldicsIAogICAgJ3JlZ2lvbl8yMGttX21pbl9lbGV2JywgJ3JlZ2lvbl8yMGttX2VsZXZfcmFuZ2UnLAogICAgJ3JlZ2lvbl81MGttX2F2Z19uZHZpJywgJ3JlZ2lvbl81MGttX2F2Z19lbGV2YXRpb24nLCAncmVnaW9uXzUwa21fYXZnX3NvaWxfbW9pc3R1cmUnLCAncmVnaW9uXzUwa21fbWF4X2VsZXYnLCAKICAgICdyZWdpb25fNTBrbV9taW5fZWxldicsICdyZWdpb25fNTBrbV9lbGV2X3JhbmdlJywKICAgICdjb3JlX3JlYWxtQWZyb3Ryb3BpYycsICdjb3JlX3JlYWxtQXVzdHJhbGFzaWEnLCAnY29yZV9yZWFsbUluZG9tYWxheWFuJywgJ2NvcmVfcmVhbG1OZWFyY3RpYycsICdjb3JlX3JlYWxtTmVvdHJvcGljJywgJ2NvcmVfcmVhbG1QYWxlYXJjdGljJykKCnR5cGVfbGFiZWxzID0gZnVuY3Rpb24ocCkgewogIGV4cGxhbmF0b3J5X2xldmVscyA9IGFsbF9leHBsYW5hdG9yaWVzW2FsbF9leHBsYW5hdG9yaWVzICVpbiUgcCRleHBsYW5hdG9yeV0KICBwJGV4cGxhbmF0b3J5IDwtIGZhY3RvcihwJGV4cGxhbmF0b3J5LCBsZXZlbHMgPSBleHBsYW5hdG9yeV9sZXZlbHMpCiAgCiAgcCR0eXBlIDwtICdSZWFsbScKICBwJHR5cGVbcCRleHBsYW5hdG9yeSAlaW4lIGMoJ2NpdHlfYXZnX25kdmknLCAnY2l0eV9hdmdfZWxldmF0aW9uJywgJ2NpdHlfYXZnX3RlbXAnLCAnY2l0eV9hdmdfbWluX21vbnRobHlfdGVtcCcsICdjaXR5X2F2Z19tYXhfbW9udGhseV90ZW1wJywgCiAgICAnY2l0eV9hdmdfbW9udGhseV90ZW1wJywgJ2NpdHlfYXZnX3JhaW5mYWxsJywgJ2NpdHlfYXZnX21heF9tb250aGx5X3JhaW5mYWxsJywgJ2NpdHlfYXZnX21pbl9tb250aGx5X3JhaW5mYWxsJywgCiAgICAnY2l0eV9hdmdfc29pbF9tb2lzdHVyZScsICdjaXR5X21heF9lbGV2JywgJ2NpdHlfbWluX2VsZXYnLCAnY2l0eV9lbGV2X3JhbmdlJyldIDwtICdDaXR5IGdlb2dyYXBoeScKICBwJHR5cGVbcCRleHBsYW5hdG9yeSAlaW4lIGMoJ3JlZ2lvbl81MGttX2F2Z19uZHZpJywgJ3JlZ2lvbl81MGttX2F2Z19lbGV2YXRpb24nLCAncmVnaW9uXzUwa21fYXZnX3NvaWxfbW9pc3R1cmUnLCAncmVnaW9uXzUwa21fbWF4X2VsZXYnLCAKICAgICdyZWdpb25fNTBrbV9taW5fZWxldicsICdyZWdpb25fNTBrbV9lbGV2X3JhbmdlJyldIDwtICdSZWdpb25hbCAoNTBrbSkgZ2VvZ3JhcGh5JwogICBwJHR5cGVbcCRleHBsYW5hdG9yeSAlaW4lIGMoJ3JlZ2lvbl8yMGttX2F2Z19uZHZpJywgJ3JlZ2lvbl8yMGttX2F2Z19lbGV2YXRpb24nLCAncmVnaW9uXzIwa21fYXZnX3NvaWxfbW9pc3R1cmUnLCAncmVnaW9uXzIwa21fbWF4X2VsZXYnLCAKICAgICdyZWdpb25fMjBrbV9taW5fZWxldicsICdyZWdpb25fMjBrbV9lbGV2X3JhbmdlJyldIDwtICdSZWdpb25hbCAoMjBrbSkgZ2VvZ3JhcGh5JwogIHAkdHlwZVtwJGV4cGxhbmF0b3J5ICVpbiUgYygnYWJzX2xhdGl0dWRlJywgJ2xhdGl0dWRlJywgJ2xvbmdpdHVkZScpXSA8LSAnU3BhdGlhbCcKICBwCn0KYGBgCgpgYGB7cn0KZXhwbGFuYXRvcnlfbGFiZWxzID0gYygKICAnaGFzX2ludHJvZHVjZWRfc3BlY2llcyc9J0hhcyBpbnRyb2R1Y2VkIHNwZWNpZXMnLCAKICAnY2l0eV9hdmdfbmR2aSc9J0F2ZXJhZ2UgTkRWSScsIAogICdjaXR5X2F2Z19lbGV2YXRpb24nPSdBdmVyYWdlIGVsZXZhdGlvbicsIAogICdjaXR5X2F2Z190ZW1wJz0nQXZlcmFnZSB0ZW1wZXJhdHVyZScsIAogICdjaXR5X2F2Z19taW5fbW9udGhseV90ZW1wJz0nQXZlcmFnZSBtaW5pbXVtIG1vbnRobHkgdGVtcGVyYXR1cmUnLCAKICAnY2l0eV9hdmdfbWF4X21vbnRobHlfdGVtcCc9J0F2ZXJhZ2UgbWF4aW11bSBtb250aGx5IHRlbXBlcmF0dXJlJywgCiAgJ2NpdHlfYXZnX21vbnRobHlfdGVtcCc9J0F2ZXJhZ2UgbW9udGhseSB0ZW1wZXJhdHVyZScsIAogICdjaXR5X2F2Z19yYWluZmFsbCc9J0F2ZXJhZ2UgcmFpbmZhbGwnLCAKICAnY2l0eV9hdmdfbWF4X21vbnRobHlfcmFpbmZhbGwnPSdBdmVyYWdlIG1heGltdW0gbW9udGhseSByYWluZmFsbCcsIAogICdjaXR5X2F2Z19taW5fbW9udGhseV9yYWluZmFsbCc9J0F2ZXJhZ2UgbWluaW11bSBtb250aGx5IHJhaW5mYWxsJywgCiAgJ2NpdHlfYXZnX3NvaWxfbW9pc3R1cmUnPSdBdmVyYWdlIHNvaWwgbW9pc3R1cmUnLCAKICAnY2l0eV9tYXhfZWxldic9J01heGltdW0gZWxldmF0aW9uJywgCiAgJ2NpdHlfbWluX2VsZXYnPSdNaW5pbXVtIGVsZXZhdGlvbicsIAogICdjaXR5X2VsZXZfcmFuZ2UnPSdFbGV2YXRpb24gcmFuZ2UnLCAKICAncmVnaW9uXzIwa21fYXZnX25kdmknPSdBdmVyYWdlIE5EVkknLCAKICAncmVnaW9uXzIwa21fYXZnX2VsZXZhdGlvbic9J0F2ZXJhZ2UgZWxldmF0aW9uJywgCiAgJ3JlZ2lvbl8yMGttX2F2Z19zb2lsX21vaXN0dXJlJz0nQXZlcmFnZSBzb2lsIG1vaXN0dXJlJywgCiAgJ3JlZ2lvbl8yMGttX21heF9lbGV2Jz0nTWF4aW11bSBlbGV2YXRpb24nLCAKICAncmVnaW9uXzIwa21fbWluX2VsZXYnPSdNaW5pbXVtIGVsZXZhdGlvbicsCiAgJ3JlZ2lvbl8yMGttX2VsZXZfcmFuZ2UnPSdFbGV2YXRpb24gcmFuZ2UnLAogICdyZWdpb25fNTBrbV9hdmdfbmR2aSc9J0F2ZXJhZ2UgTkRWSScsCiAgJ3JlZ2lvbl81MGttX2F2Z19lbGV2YXRpb24nPSdBdmVyYWdlIGVsZXZhdGlvbicsCiAgJ3JlZ2lvbl81MGttX2F2Z19zb2lsX21vaXN0dXJlJz0nQXZlcmFnZSBzb2lsIG1vaXN0dXJlJywgCiAgJ3JlZ2lvbl81MGttX21heF9lbGV2Jz0nTWF4aW11bSBlbGV2YXRpb24nLAogICdyZWdpb25fNTBrbV9taW5fZWxldic9J01pbmltdW0gZWxldmF0aW9uJywgCiAgJ3JlZ2lvbl81MGttX2VsZXZfcmFuZ2UnPSdFbGV2YXRpb24gcmFuZ2UnLAogICdhYnNfbGF0aXR1ZGUnID0gJ0Fic29sdXRlIGxhdGl0dWRlJywKICAnbGF0aXR1ZGUnID0gJ0xhdGl0dWRlJywKICAnbG9uZ2l0dWRlJyA9ICdMb25naXR1ZGUnLAogICdjb3JlX3JlYWxtQWZyb3Ryb3BpYycgPSAnQWZyb3Ryb3BpY2FsJywgCiAgJ2NvcmVfcmVhbG1BdXN0cmFsYXNpYScgPSAnQXVzdGFsaWFzaWFuJywgCiAgJ2NvcmVfcmVhbG1JbmRvbWFsYXlhbicgPSAnSW5kb21hbGF5YW4nLCAKICAnY29yZV9yZWFsbU5lYXJjdGljJyA9ICdOZWFyY3RpYycsIAogICdjb3JlX3JlYWxtTmVvdHJvcGljJyA9ICdOZW90cm9waWNhbCcsCiAgJ2NvcmVfcmVhbG1QYWxlYXJjdGljJyA9ICdQYWxlYXJjdGljJywKICAnY29yZV9yZWFsbU9jZWFuaWEnID0gJ09jZWFuaWNhbCcpCmBgYAoKIyMgSGVscGVyIHBsb3QgZnVuY3Rpb25zCmBgYHtyfQpnZW9tX25vcm1hbGlzZWRfaGlzdG9ncmFtID0gZnVuY3Rpb24obmFtZSwgZ2csIGxlZ2VuZC5wb3NpdGlvbiA9ICJib3R0b20iKSB7CiAgZ2cgKyAKICAgIGdlb21faGlzdG9ncmFtKGFlcyhmaWxsID0gY29yZV9yZWFsbSksIGJpbndpZHRoID0gMC4xLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsKICAgIGdlb21fdmxpbmUoYWVzKHhpbnRlcmNlcHQgPSAwLjUpLCBjb2xvciA9ICIjMDAwMDAwIiwgc2l6ZSA9IDAuNCkgKwogICAgZ2VvbV92bGluZShhZXMoeGludGVyY2VwdCA9IDApLCBjb2xvciA9ICIjMDAwMDAwIiwgc2l6ZSA9IDAuMiwgbGluZXR5cGUgPSAiZGFzaGVkIikgKwogICAgZ2VvbV92bGluZShhZXMoeGludGVyY2VwdCA9IDEpLCBjb2xvciA9ICIjMDAwMDAwIiwgc2l6ZSA9IDAuMiwgbGluZXR5cGUgPSAiZGFzaGVkIikgKyAKICAgIHlsYWIoIk51bWJlciBvZiBjaXRpZXMiKSArIHhsYWIoIk5vcm1hbGlzZWQgUmVzcG9uc2UiKSArIHlsaW0oYygwLCA3MCkpICsKICAgIGxhYnModGl0bGUgPSBuYW1lLCBmaWxsID0gJ1JlYWxtJykgKwogICAgdGhlbWVfYncoKSArCiAgICB0aGVtZShsZWdlbmQucG9zaXRpb249bGVnZW5kLnBvc2l0aW9uKQp9CmBgYAoKYGBge3J9Cmdlb21fbWFwID0gZnVuY3Rpb24obWFwX3NmLCB0aXRsZSkgewogIG5vcm1fbW50ZF9hbmFseXNpc19nZW8gPSBnZ3Bsb3QoKSArIAogICAgZ2VvbV9zZihkYXRhID0gd29ybGRfbWFwLCBhZXMoZ2VvbWV0cnkgPSBnZW9tZXRyeSkpICsKICAgIG1hcF9zZiArCiAgICBub3JtYWxpc2VkX2NvbG91cnNfc2NhbGUgKwogICAgbGFicyhjb2xvdXIgPSAnTm9ybWFsaXNlZFxuUmVzcG9uc2UnKSArCiAgICB0aGVtZV9idygpICsKICAgIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0iYm90dG9tIikKfQpgYGAKCiMjIEhlbHBlciBEcmVkZ2UgZnVuY3Rpb25zCmBgYHtyfQojIFRha2VuIGZyb20gTXVNSU4gcGFja2FnZQojIGh0dHBzOi8vcmRyci5pby9jcmFuL011TUluL3NyYy9SL2F2ZXJhZ2luZy5SCiMgaHR0cHM6Ly9yZHJyLmlvL2NyYW4vTXVNSW4vc3JjL1IvbW9kZWwuYXZnLlIKCi5jb2VmYXJyLmF2ZyA8LQogIGZ1bmN0aW9uKGNmYXJyLCB3ZWlnaHQsIHJldmlzZWQudmFyLCBmdWxsLCBhbHBoYSkgewkKICAgIHdlaWdodCA8LSB3ZWlnaHQgLyBzdW0od2VpZ2h0KQogICAgbkNvZWYgPC0gZGltKGNmYXJyKVszTF0KICAgIGlmKGZ1bGwpIHsKICAgICAgbmFzIDwtIGlzLm5hKGNmYXJyWywgMUwsIF0pICYgaXMubmEoY2ZhcnJbLCAyTCwgXSkKICAgICAgY2ZhcnJbLCAxTCwgXVtuYXNdIDwtIGNmYXJyWywgMkwsIF1bbmFzXSA8LSAwCiAgICAgICNjZmFyclssIDFMOjJMLCBdW2lzLm5hKGNmYXJyWywgMUw6MkwsIF0pXSA8LSAwCiAgICAgIGlmKCFhbGwoaXMubmEoY2ZhcnJbLCAzTCwgXSkpKQogICAgICAgIGNmYXJyWyAsM0wsIF1baXMubmEoY2ZhcnJbICwgM0wsIF0pXSA8LSBJbmYKICAgIH0KICAgIAogICAgYXZnY29lZiA8LSBhcnJheShkaW0gPSBjKG5Db2VmLCA1TCksCiAgICAgICAgICAgICAgICAgICAgIGRpbW5hbWVzID0gbGlzdChkaW1uYW1lcyhjZmFycilbWzNMXV0sIGMoIkVzdGltYXRlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU3RkLiBFcnJvciIsICJBZGp1c3RlZCBTRSIsICJMb3dlciBDSSIsICJVcHBlciBDSSIpKSkKICAgIGZvcihpIGluIHNlcV9sZW4obkNvZWYpKQogICAgICBhdmdjb2VmW2ksIF0gPC0gcGFyLmF2ZyhjZmFyclssIDFMLCBpXSwgY2ZhcnJbLCAyTCwgaV0sIHdlaWdodCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGYgPSBjZmFyclssIDNMLCBpXSwgYWxwaGEgPSBhbHBoYSwgcmV2aXNlZC52YXIgPSByZXZpc2VkLnZhcikKICAgIAogICAgYXZnY29lZltpcy5uYW4oYXZnY29lZildIDwtIE5BCiAgICByZXR1cm4oYXZnY29lZikKICB9CgoubWFrZWNvZWZtYXQgPC0gZnVuY3Rpb24oY2YpIHsKICBuby5hc2UgPC0gYWxsKGlzLm5hKGNmWywgM0xdKSkKICB6IDwtIGFicyhjZlssIDFMXSAvIGNmWywgaWYobm8uYXNlKSAyTCBlbHNlIDNMXSkKICBwdmFsIDwtIDIgKiBwbm9ybSh6LCBsb3dlci50YWlsID0gRkFMU0UpCiAgY2JpbmQoY2ZbLCBpZihuby5hc2UpIDFMOjJMIGVsc2UgMUw6M0wsIGRyb3AgPSBGQUxTRV0sCiAgICAgICAgYHogdmFsdWVgID0geiwgYFByKD58enwpYCA9IHphcHNtYWxsKHB2YWwpKQp9CgojIEdlbmVyYXRlIG1vZGVsIHNlbGVjdGlvbnMgdXNpbmcgbG1lciwgZHJlZGdlLCBhbmQgbW9kZWwuYXZnCiMgYGZvcnVtbGFgIDogYSB0d28tc2lkZWQgbGluZWFyIGZvcm11bGEgb2JqZWN0IGRlc2NyaWJpbmcgYm90aCB0aGUgZml4ZWQtZWZmZWN0cyBhbmQgcmFuZG9tLWVmZmVjdHMgcGFydCBvZiB0aGUgbW9kZWwKIyBgZGF0YWAgOiB0aGUgZGF0YSBmcmFtZSBjb250YWluaW5nIHRoZSB2YXJpYWJsZXMgZnJvbSB0aGUgZm9ybXVsYQojIGBhaWNfZGVsdGFgIDogdGhlIEFJQyBkZWx0YSB0byB1c2UgZm9yIHNlbGVjdGluZyBtb2RlbHMgaW4gbW9kZWwgYXZlcmFnZQptb2RlbF9hdmVyYWdlIDwtIGZ1bmN0aW9uKGZvcm11bGEsIGRhdGEsIGFpY19kZWx0YSA9IDIwKSB7CiAgbW9kZWwgPC0gbG0oCiAgICBmb3JtdWxhLAogICAgZGF0YT1kYXRhCiAgKQogIGRyZWRnZV9yZXN1bHQgPC0gZHJlZGdlKG1vZGVsKQogIHN1bW1hcnkobW9kZWwuYXZnKGRyZWRnZV9yZXN1bHQsIHN1YnNldCA9IGRlbHRhIDwgYWljX2RlbHRhKSkKfQoKIyBDcmVhdGUgYSBzdW1tYXJ5IGRhdGEgZnJhbWUgY29udGFpbmluZyB0aGUgc2VsZWN0ZWQgdmFyaWFibGVzIGZyb20gYSBtb2RlbAojIGBtb2RlbF9zdW1gIDogVGhlIG1vZGVsIHN1bW1hcnkgb3V0cHV0IGZyb20gYG1vZGVsX2F2ZXJhZ2VgCm1vZGVsX3N1bW1hcnkgPC0gZnVuY3Rpb24obW9kZWxfc3VtKSB7CiAgLmNvbHVtbl9uYW1lIDwtIGZ1bmN0aW9uKHBvc3RmaXgpIHsKICAgIHBvc3RmaXgKICB9CiAgCiAgIyBqdXN0IHJldHVybiB0aGUgZXN0aW1hdGUgYW5kIHAgdmFsdWUKICB3ZWlnaHQgPC0gbW9kZWxfc3VtJG1zVGFibGVbLCA1TF0KICAKICBjb2VmbWF0LmZ1bGwgPC0gYXMuZGF0YS5mcmFtZSgubWFrZWNvZWZtYXQoLmNvZWZhcnIuYXZnKG1vZGVsX3N1bSRjb2VmQXJyYXksIHdlaWdodCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGF0dHIobW9kZWxfc3VtLCAicmV2aXNlZC52YXIiKSwgVFJVRSwgMC4wNSkpKQogIAogIGNvZWZtYXQuc3Vic2V0IDwtCiAgICBhcy5kYXRhLmZyYW1lKC5tYWtlY29lZm1hdCguY29lZmFyci5hdmcobW9kZWxfc3VtJGNvZWZBcnJheSwgd2VpZ2h0LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGF0dHIobW9kZWxfc3VtLCAicmV2aXNlZC52YXIiKSwgRkFMU0UsIDAuMDUpKSkKICAKICAKICBjb2VmbWF0LnN1YnNldCA8LSBjb2VmbWF0LnN1YnNldFstYygxKSwgYygxLCAyLCA1KV0KICBuYW1lcyhjb2VmbWF0LnN1YnNldCkgPC0gYyguY29sdW1uX25hbWUoImVzdGltYXRlIiksIC5jb2x1bW5fbmFtZSgiZXJyb3IiKSwgLmNvbHVtbl9uYW1lKCJwIikpCiAgY29lZm1hdC5zdWJzZXQgPC0gdGliYmxlOjpyb3duYW1lc190b19jb2x1bW4oY29lZm1hdC5zdWJzZXQsICJleHBsYW5hdG9yeSIpCiAgY29lZm1hdC5zdWJzZXQkbW9kZWwgPSAnc3Vic2V0JwogIAogIGNvZWZtYXQuZnVsbCA8LSBjb2VmbWF0LmZ1bGxbLWMoMSksIGMoMSwgMiwgNSldCiAgbmFtZXMoY29lZm1hdC5mdWxsKSA8LSBjKC5jb2x1bW5fbmFtZSgiZXN0aW1hdGUiKSwgLmNvbHVtbl9uYW1lKCJlcnJvciIpLCAuY29sdW1uX25hbWUoInAiKSkKICBjb2VmbWF0LmZ1bGwgPC0gdGliYmxlOjpyb3duYW1lc190b19jb2x1bW4oY29lZm1hdC5mdWxsLCAiZXhwbGFuYXRvcnkiKQogIGNvZWZtYXQuZnVsbCRtb2RlbCA9ICdmdWxsJwogIAogIHJiaW5kKGNvZWZtYXQuZnVsbCwgY29lZm1hdC5zdWJzZXQpCn0KYGBgCgpgYGB7cn0KZm9ybXVsYV9mcm9tX3ZzdXJwID0gZnVuY3Rpb24ocHJlZGljdG9ycywgZGVwZW5kZW50LCB2c3VycF9yZXN1bHQpIHsKICBhcy5mb3JtdWxhKHBhc3RlKGRlcGVuZGVudCwgcGFzdGUobmFtZXMocHJlZGljdG9yc1ssdnN1cnBfcmVzdWx0JHZhcnNlbGVjdC5pbnRlcnBdKSwgY29sbGFwc2U9IisiKSwgc2VwID0gIn4iKSkKfQpgYGAKCmBgYHtyfQpwbG90X3ZzdXJwX3Jlc3VsdCA9IGZ1bmN0aW9uKHJlc3VsdF90YWJsZSkgewogIHBsb3QgPSByZXN1bHRfdGFibGVbcmVzdWx0X3RhYmxlJG1vZGVsID09ICdmdWxsJyxdCiAgcGxvdCA9IHR5cGVfbGFiZWxzKHBsb3QpCgogIGdncGxvdChwbG90LCBhZXMoeT1leHBsYW5hdG9yeSwgeD1lc3RpbWF0ZSwgY29sb3VyID0gdHlwZSkpICsgCiAgICBnZW9tX2xpbmUoKSArCiAgICBnZW9tX3BvaW50KCkrCiAgICBnZW9tX2Vycm9yYmFyKGFlcyh4bWluPWVzdGltYXRlLWVycm9yLCB4bWF4PWVzdGltYXRlK2Vycm9yKSwgd2lkdGg9LjIsCiAgICAgICAgICAgICAgICAgICBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSgwLjA1KSkgKwogICAgc2NhbGVfeV9kaXNjcmV0ZSgKICAgICAgbGltaXRzID0gcmV2KGxldmVscyhwbG90JGV4cGxhbmF0b3J5KSksIAogICAgICBsYWJlbHMgPSBleHBsYW5hdG9yeV9sYWJlbHMpICsKICAgIHRoZW1lX2J3KCkgKwogICAgZ2VvbV92bGluZSh4aW50ZXJjZXB0PTAsIGxpbmV0eXBlPSJkb3R0ZWQiKSArCiAgICBndWlkZXMoY29sb3VyPWd1aWRlX2xlZ2VuZCh0aXRsZT0iUHJlZGljdG9yIHR5cGUiKSkgKyB4bGFiKCdJbmNyZWFzZSBmcm9tIDAgKGhhYml0YXQgZmlsdGVyaW5nKVxudG8gMSAoY29tcGV0aXRpdmUgZXhjbHVzaW9uKVxuwrEgU3RhbmRhcmQgRXJyb3InKSArIHlsYWIoJ1ByZWRpY3RvcicpICsKICAgIHRoZW1lKGxlZ2VuZC5qdXN0aWZpY2F0aW9uID0gInRvcCIpCn0KYGBgCgojIyBNTlRECmBgYHtyfQpub3JtX21udGRfYW5hbHlzaXNfcGxvdCA9IGdlb21fbm9ybWFsaXNlZF9oaXN0b2dyYW0oCiAgJ01OVEQnLCAKICBnZ3Bsb3QoYW5hbHlzaXNfZGF0YSwgYWVzKG1udGRfbm9ybWFsaXNlZCkpCikKbm9ybV9tbnRkX2FuYWx5c2lzX3Bsb3QKYGBgCgpgYGB7cn0Kbm9ybV9tbnRkX2FuYWx5c2lzX2dlb19wbG90ID0gZ2VvbV9tYXAoZ2VvbV9zZihkYXRhID0gYW5hbHlzaXNfZGF0YSwgYWVzKGNvbG9yID0gbW50ZF9ub3JtYWxpc2VkLCBnZW9tZXRyeSA9IGdlb21ldHJ5KSksICdNTlREJykKbm9ybV9tbnRkX2FuYWx5c2lzX2dlb19wbG90CmBgYAoKYGBge3J9Cm5vcm1fbW50ZF9hbmFseXNpc19kYXRhID0gbW9kZWxfZGF0YShhbmFseXNpc19kYXRhWyFpcy5uYShhbmFseXNpc19kYXRhJG1udGRfbm9ybWFsaXNlZCksXSwgJ21udGRfbm9ybWFsaXNlZCcpCm5vcm1fbW50ZF9hbmFseXNpc19wcmVkaWN0b3JzID0gbm9ybV9tbnRkX2FuYWx5c2lzX2RhdGFbLC0xXQpub3JtX21udGRfYW5hbHlzaXNfaW50ZXJwID0gVlNVUkYoeCA9IG5vcm1fbW50ZF9hbmFseXNpc19wcmVkaWN0b3JzLCB5ID0gbm9ybV9tbnRkX2FuYWx5c2lzX2RhdGEkbW50ZF9ub3JtYWxpc2VkKQpuYW1lcyhub3JtX21udGRfYW5hbHlzaXNfcHJlZGljdG9yc1ssbm9ybV9tbnRkX2FuYWx5c2lzX2ludGVycCR2YXJzZWxlY3QuaW50ZXJwXSkKYGBgCgpgYGB7cn0Kbm9ybV9tbnRkX2FuYWx5c2lzX2Zvcm11bGEgPSBmb3JtdWxhX2Zyb21fdnN1cnAobm9ybV9tbnRkX2FuYWx5c2lzX3ByZWRpY3RvcnMsICJtbnRkX25vcm1hbGlzZWQiLCBub3JtX21udGRfYW5hbHlzaXNfaW50ZXJwKQpub3JtX21udGRfYW5hbHlzaXNfcmVzdWx0IDwtIG1vZGVsX2F2ZXJhZ2Uobm9ybV9tbnRkX2FuYWx5c2lzX2Zvcm11bGEsIG5vcm1fbW50ZF9hbmFseXNpc19kYXRhKQpub3JtX21udGRfYW5hbHlzaXNfcmVzdWx0X3RhYmxlID0gbW9kZWxfc3VtbWFyeShub3JtX21udGRfYW5hbHlzaXNfcmVzdWx0KQpub3JtX21udGRfYW5hbHlzaXNfcmVzdWx0X3RhYmxlCmBgYAoKYGBge3J9Cm5vcm1fbW50ZF9hbmFseXNpc19wcmVkX3Bsb3QgPSBwbG90X3ZzdXJwX3Jlc3VsdChub3JtX21udGRfYW5hbHlzaXNfcmVzdWx0X3RhYmxlKQpub3JtX21udGRfYW5hbHlzaXNfcHJlZF9wbG90CmBgYAoKIyMgR2FwZSB3aWR0aCAtIEZEaXYKYGBge3J9Cm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3Bsb3QgPSBnZW9tX25vcm1hbGlzZWRfaGlzdG9ncmFtKAogICdHYXBlIEZEaXYnLCAKICBnZ3Bsb3QoYW5hbHlzaXNfZGF0YSwgYWVzKGdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkKSkKKQpub3JtX2dhcGVfZmRpdl9hbmFseXNpc19wbG90CmBgYAoKYGBge3J9Cm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX2dlb19wbG90ID0gZ2VvbV9tYXAoZ2VvbV9zZihkYXRhID0gYW5hbHlzaXNfZGF0YSwgYWVzKGNvbG9yID0gZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQsIGdlb21ldHJ5ID0gZ2VvbWV0cnkpKSwgJ0dhcGUgV2lkdGggRkRpdicpCm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX2dlb19wbG90CmBgYAoKYGBge3J9Cm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX2RhdGEgPSBtb2RlbF9kYXRhKGFuYWx5c2lzX2RhdGFbIWlzLm5hKGFuYWx5c2lzX2RhdGEkZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQpLF0sICdnYXBlX3dpZHRoX2ZkaXZfbm9ybWFsaXNlZCcpCm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnMgPSBub3JtX2dhcGVfZmRpdl9hbmFseXNpc19kYXRhWywtMV0Kbm9ybV9nYXBlX2ZkaXZfYW5hbHlzaXNfaW50ZXJwID0gVlNVUkYoeCA9IG5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnMsIHkgPSBub3JtX2dhcGVfZmRpdl9hbmFseXNpc19kYXRhJGdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkKQpuYW1lcyhub3JtX2dhcGVfZmRpdl9hbmFseXNpc19wcmVkaWN0b3JzWyxub3JtX2dhcGVfZmRpdl9hbmFseXNpc19pbnRlcnAkdmFyc2VsZWN0LmludGVycF0pCmBgYAoKYGBge3J9Cm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX2Zvcm11bGEgPSBmb3JtdWxhX2Zyb21fdnN1cnAobm9ybV9nYXBlX2ZkaXZfYW5hbHlzaXNfcHJlZGljdG9ycywgImdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkIiwgbm9ybV9nYXBlX2ZkaXZfYW5hbHlzaXNfaW50ZXJwKQpub3JtX2dhcGVfZmRpdl9hbmFseXNpc19yZXN1bHQgPC0gbW9kZWxfYXZlcmFnZShub3JtX2dhcGVfZmRpdl9hbmFseXNpc19mb3JtdWxhLCBub3JtX2dhcGVfZmRpdl9hbmFseXNpc19kYXRhKQpub3JtX2dhcGVfZmRpdl9hbmFseXNpc19yZXN1bHRfdGFibGUgPSBtb2RlbF9zdW1tYXJ5KG5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3Jlc3VsdCkKbm9ybV9nYXBlX2ZkaXZfYW5hbHlzaXNfcmVzdWx0X3RhYmxlCmBgYAoKYGBge3J9Cm5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3ByZWRfcGxvdCA9IHBsb3RfdnN1cnBfcmVzdWx0KG5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3Jlc3VsdF90YWJsZSkKbm9ybV9nYXBlX2ZkaXZfYW5hbHlzaXNfcHJlZF9wbG90CmBgYAoKCiMjIEhXSSAtIEZEaXYKYGBge3J9Cm5vcm1faHdpX2ZkaXZfbG9jb19hbmFseXNpc19wbG90ID0gZ2VvbV9ub3JtYWxpc2VkX2hpc3RvZ3JhbSgKICAnSFdJIEZEaXYnLCAKICBnZ3Bsb3QoYW5hbHlzaXNfZGF0YSwgYWVzKGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCkpCikKbm9ybV9od2lfZmRpdl9sb2NvX2FuYWx5c2lzX3Bsb3QKYGBgCgpgYGB7cn0Kbm9ybV9od2lfZmRpdl9hbmFseXNpc19nZW9fcGxvdCA9IGdlb21fbWFwKGdlb21fc2YoZGF0YSA9IGFuYWx5c2lzX2RhdGEsIGFlcyhjb2xvciA9IGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCwgZ2VvbWV0cnkgPSBnZW9tZXRyeSkpLCAnSFdJIEZEaXYnKQpub3JtX2h3aV9mZGl2X2FuYWx5c2lzX2dlb19wbG90CmBgYAoKYGBge3J9Cm5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfZGF0YSA9IG1vZGVsX2RhdGEoYW5hbHlzaXNfZGF0YVshaXMubmEoYW5hbHlzaXNfZGF0YSRoYW5kd2luZ19pbmRleF9mZGl2X25vcm1hbGlzZWQpLF0sICdoYW5kd2luZ19pbmRleF9mZGl2X25vcm1hbGlzZWQnKQpub3JtX2h3aV9mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnMgPSBub3JtX2h3aV9mZGl2X2FuYWx5c2lzX2RhdGFbLC0xXQpub3JtX2h3aV9mZGl2X2FuYWx5c2lzX2ludGVycCA9IFZTVVJGKHggPSBub3JtX2h3aV9mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnMsIHkgPSBub3JtX2h3aV9mZGl2X2FuYWx5c2lzX2RhdGEkaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkKQpuYW1lcyhub3JtX2h3aV9mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnNbLG5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfaW50ZXJwJHZhcnNlbGVjdC5pbnRlcnBdKQpgYGAKCgpgYGB7cn0Kbm9ybV9od2lfZmRpdl9hbmFseXNpc19mb3JtdWxhID0gZm9ybXVsYV9mcm9tX3ZzdXJwKG5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfcHJlZGljdG9ycywgImhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCIsIG5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfaW50ZXJwKQpub3JtX2h3aV9mZGl2X2FuYWx5c2lzX3Jlc3VsdCA8LSBtb2RlbF9hdmVyYWdlKG5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfZm9ybXVsYSwgbm9ybV9od2lfZmRpdl9hbmFseXNpc19kYXRhKQpub3JtX2h3aV9mZGl2X2FuYWx5c2lzX3Jlc3VsdF90YWJsZSA9IG1vZGVsX3N1bW1hcnkobm9ybV9od2lfZmRpdl9hbmFseXNpc19yZXN1bHQpCm5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfcmVzdWx0X3RhYmxlCmBgYAoKYGBge3J9Cm5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfcHJlZF9wbG90ID0gcGxvdF92c3VycF9yZXN1bHQobm9ybV9od2lfZmRpdl9hbmFseXNpc19yZXN1bHRfdGFibGUpCm5vcm1faHdpX2ZkaXZfYW5hbHlzaXNfcHJlZF9wbG90CmBgYAoKCiMjIE1hc3MgLSBGRGl2CmBgYHtyfQpub3JtX21hc3NfZmRpdl9sb2NvX2FuYWx5c2lzX3Bsb3QgPSBnZW9tX25vcm1hbGlzZWRfaGlzdG9ncmFtKAogICdNYXNzIEZEaXYnLCAKICBnZ3Bsb3QoYW5hbHlzaXNfZGF0YSwgYWVzKG1hc3NfZmRpdl9ub3JtYWxpc2VkKSkKKQpub3JtX21hc3NfZmRpdl9sb2NvX2FuYWx5c2lzX3Bsb3QKYGBgCgpgYGB7cn0Kbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfZ2VvX3Bsb3QgPSBnZW9tX21hcChnZW9tX3NmKGRhdGEgPSBhbmFseXNpc19kYXRhLCBhZXMoY29sb3IgPSBtYXNzX2ZkaXZfbm9ybWFsaXNlZCwgZ2VvbWV0cnkgPSBnZW9tZXRyeSkpLCAnTWFzcyBGRGl2JykKbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfZ2VvX3Bsb3QKYGBgCgpgYGB7cn0Kbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfZGF0YSA9IG1vZGVsX2RhdGEoYW5hbHlzaXNfZGF0YVshaXMubmEoYW5hbHlzaXNfZGF0YSRtYXNzX2ZkaXZfbm9ybWFsaXNlZCksXSwgJ21hc3NfZmRpdl9ub3JtYWxpc2VkJykKbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfcHJlZGljdG9ycyA9IG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2RhdGFbLC0xXQpub3JtX21hc3NfZmRpdl9hbmFseXNpc19pbnRlcnAgPSBWU1VSRih4ID0gbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfcHJlZGljdG9ycywgeSA9IG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2RhdGEkbWFzc19mZGl2X25vcm1hbGlzZWQpCm5hbWVzKG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX3ByZWRpY3RvcnNbLG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2ludGVycCR2YXJzZWxlY3QuaW50ZXJwXSkKYGBgCgpgYGB7cn0Kbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfZm9ybXVsYSA9IGZvcm11bGFfZnJvbV92c3VycChub3JtX21hc3NfZmRpdl9hbmFseXNpc19wcmVkaWN0b3JzLCAibWFzc19mZGl2X25vcm1hbGlzZWQiLCBub3JtX21hc3NfZmRpdl9hbmFseXNpc19pbnRlcnApCm5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX3Jlc3VsdCA8LSBtb2RlbF9hdmVyYWdlKG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2Zvcm11bGEsIG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2RhdGEpCm5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX3Jlc3VsdF90YWJsZSA9IG1vZGVsX3N1bW1hcnkobm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfcmVzdWx0KQpub3JtX21hc3NfZmRpdl9hbmFseXNpc19yZXN1bHRfdGFibGUKYGBgCgpgYGB7cn0Kbm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfcHJlZF9wbG90ID0gcGxvdF92c3VycF9yZXN1bHQobm9ybV9tYXNzX2ZkaXZfYW5hbHlzaXNfcmVzdWx0X3RhYmxlKQpub3JtX21hc3NfZmRpdl9hbmFseXNpc19wcmVkX3Bsb3QKYGBgCgojIENyZWF0ZSBwbG90IG9mIGRpZmZlcmVuY2VzIGluIHByb2Nlc3MgcmVzcG9uc2UKYGBge3J9CnByZWRfbGVnZW5kIDwtIGdldF9sZWdlbmQoCiAgIyBjcmVhdGUgc29tZSBzcGFjZSB0byB0aGUgbGVmdCBvZiB0aGUgbGVnZW5kCiAgbm9ybV9od2lfZmRpdl9hbmFseXNpc19wcmVkX3Bsb3QgKyB0aGVtZShsZWdlbmQuYm94Lm1hcmdpbiA9IG1hcmdpbigwLCAwLCAwLCAxMikpCikKZ2VvX2xlZ2VuZCA8LSBnZXRfbGVnZW5kKAogICMgY3JlYXRlIHNvbWUgc3BhY2UgdG8gdGhlIGxlZnQgb2YgdGhlIGxlZ2VuZAogIG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX2dlb19wbG90ICsgdGhlbWUobGVnZW5kLmJveC5tYXJnaW4gPSBtYXJnaW4oMCwgMCwgMCwgMTIpKQopCgpwbG90X2dyaWQoCiAgcGxvdF9ncmlkKAogICAgbm9ybV9tbnRkX2FuYWx5c2lzX2dlb19wbG90ICsgdGhlbWUobGVnZW5kLnBvc2l0aW9uPSJub25lIiksIAogICAgbm9ybV9tbnRkX2FuYWx5c2lzX3ByZWRfcGxvdCArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpLCAKICAgIG5yb3cgPSAxCiAgKSArIGRyYXdfbGFiZWwoIk1OVEQiLCBzaXplID0gMTYsIGFuZ2xlID0gOTAsIHggPSAwLjAxLCB5ID0gMC41KSwKICBwbG90X2dyaWQoCiAgICBub3JtX2dhcGVfZmRpdl9hbmFseXNpc19nZW9fcGxvdCArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpLCAKICAgIG5vcm1fZ2FwZV9mZGl2X2FuYWx5c2lzX3ByZWRfcGxvdCArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpLCAKICAgIG5yb3cgPSAxCiAgKSArIGRyYXdfbGFiZWwoIkdhcGUiLCBzaXplID0gMTYsIGFuZ2xlID0gOTAsIHggPSAwLjAxLCB5ID0gMC41KSwKICBwbG90X2dyaWQoCiAgICBub3JtX2h3aV9mZGl2X2FuYWx5c2lzX2dlb19wbG90ICsgdGhlbWUobGVnZW5kLnBvc2l0aW9uPSJub25lIiksIAogICAgbm9ybV9od2lfZmRpdl9hbmFseXNpc19wcmVkX3Bsb3QgKyB0aGVtZShsZWdlbmQucG9zaXRpb249Im5vbmUiKSwgCiAgICBucm93ID0gMQogICkgKyBkcmF3X2xhYmVsKCJIV0kiLCBzaXplID0gMTYsIGFuZ2xlID0gOTAsIHggPSAwLjAxLCB5ID0gMC41KSwKICBwbG90X2dyaWQoCiAgICBub3JtX21hc3NfZmRpdl9hbmFseXNpc19nZW9fcGxvdCArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpLCAKICAgIG5vcm1fbWFzc19mZGl2X2FuYWx5c2lzX3ByZWRfcGxvdCArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpLCAKICAgIG5yb3cgPSAxCiAgKSArIGRyYXdfbGFiZWwoIk1hc3MiLCBzaXplID0gMTYsIGFuZ2xlID0gOTAsIHggPSAwLjAxLCB5ID0gMC41KSwgCiAgcGxvdF9ncmlkKAogICAgZ2VvX2xlZ2VuZCwgCiAgICBwcmVkX2xlZ2VuZCwgCiAgICBucm93ID0gMQogICksCiAgbnJvdyA9IDUKKQpnZ3NhdmUoZmlsZW5hbWUoRklHVVJFU19PVVRQVVRfRElSLCAncHJvY2Vzc19yZXNwb25zZS5qcGcnKSwgd2lkdGggPSAyNTAwLCBoZWlnaHQ9NTAwMCwgdW5pdHMgPSAncHgnKQpgYGAKCgojIENvbXBhcmUgbWV0cmljcyBhZ2FpbnN0IGVhY2ggb3RoZXIKYGBge3J9CmdncGxvdChhbmFseXNpc19kYXRhLCBhZXMoeCA9IGdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkLCB5ID0gbW50ZF9ub3JtYWxpc2VkLCBjb2xvdXIgPSBjb3JlX3JlYWxtKSkgKyAKICBnZW9tX3BvaW50KCkgKwogIHlsYWIoIk1OVEQiKSArIAogIHhsYWIoIkdhcGUgV2lkdGggRkRpdiIpICsKICB0aGVtZV9idygpICsgbGFicyhjb2xvciA9ICJSZWFsbSIpCmBgYAoKYGBge3J9CmdncGxvdChhbmFseXNpc19kYXRhLCBhZXMoeCA9IGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCwgeSA9IG1udGRfbm9ybWFsaXNlZCwgY29sb3VyID0gY29yZV9yZWFsbSkpICsgCiAgZ2VvbV9wb2ludCgpICsKICB5bGFiKCJNTlREIikgKyAKICB4bGFiKCJIV0kgRkRpdiIpICsKICB0aGVtZV9idygpICsgbGFicyhjb2xvciA9ICJSZWFsbSIpCmBgYAoKYGBge3J9CmdncGxvdChhbmFseXNpc19kYXRhLCBhZXMoeCA9IGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCwgeSA9IGdhcGVfd2lkdGhfZmRpdl9ub3JtYWxpc2VkLCBjb2xvdXIgPSBjb3JlX3JlYWxtKSkgKyAKICBnZW9tX3BvaW50KCkgKwogIHlsYWIoIkdhcGUgV2lkdGggRkRpdiIpICsgCiAgeGxhYigiSFdJIEZEaXYiKSArCiAgdGhlbWVfYncoKSArIGxhYnMoY29sb3IgPSAiUmVhbG0iKQpgYGAKCmBgYHtyfQptbnRkX2ZkaXZfYW5hbHlzaXMgPSBhbmFseXNpc19kYXRhICU+JSAKICBkcGx5cjo6c2VsZWN0KGNpdHlfaWQsICBtbnRkX25vcm1hbGlzZWQsIGhhbmR3aW5nX2luZGV4X2ZkaXZfbm9ybWFsaXNlZCwgZ2FwZV93aWR0aF9mZGl2X25vcm1hbGlzZWQpICU+JQogIGxlZnRfam9pbihjb21tdW5pdHlfc3VtbWFyeSkgJT4lCiAgbXV0YXRlKHVyYmFuX3Bvb2xfcGVyYyA9IHVyYmFuX3Bvb2xfc2l6ZSAqIDEwMCAvIHJlZ2lvbmFsX3Bvb2xfc2l6ZSkKbW50ZF9mZGl2X2FuYWx5c2lzCmBgYAoKYGBge3J9CmdncGFpcnMobW50ZF9mZGl2X2FuYWx5c2lzICU+JSBkcGx5cjo6c2VsZWN0KG1udGRfbm9ybWFsaXNlZCwgaGFuZHdpbmdfaW5kZXhfZmRpdl9ub3JtYWxpc2VkLCBnYXBlX3dpZHRoX2ZkaXZfbm9ybWFsaXNlZCwgcmVnaW9uYWxfcG9vbF9zaXplLCB1cmJhbl9wb29sX3NpemUsIHVyYmFuX3Bvb2xfcGVyYykpCmdnc2F2ZShmaWxlbmFtZShGSUdVUkVTX09VVFBVVF9ESVIsICdhcHBlbmRpeF9ub3JtYWxpc2VkX2NvcnJlbGF0aW9uLmpwZycpKQpgYGAKCgo=